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Abstract 

For a field k with an automorphism a and a derivation 5, we intro- 
duce the notion of houvilhan sohitions of hnear difference-differential 
systems {cr(i^) = AY,S{Y) = BY} over k and characterize the exis- 
tence of liouvillian solutions in terms of the Galois group of the systems. 
We will give an algorithm to decide whether such a system has liou- 
villian solutions when k ~ C(a;, t\ a{x) = .t + 1, (5 = ^ and the size of 
the system is a prime. 

Contents 

1 Introduction 

2 Galois Theory [4 

2.1 Picard-Vessiot extensions and Galois groups J4 

2.2 Liouvillian solutions 

2.3 From e^R to R __ 

2.4 From C(x,t) to C(t)(j;) 



*Key Laboratory of Mathematics Mechanization, Institute of Systems Science, AMSS, 
CAS, Beijing 100190, China, rfeng2@ncsu.edu. The author is supported by NKBRPC 
2004CB318000 and NSFC10671200. The worft was done during a stay of the first autlior 
at Department of Mathematics, North Carolina State University (NCSU). The hospitality 
at NCSU is gratefully acknowledged. 

^North Carolina State University, Department of Mathematics, Box 8205, Raleigh, 
North Carolina 27695-8205, USA, singer@math.ncsu.edu. This material is based upon 
worli supported by the National Science Foundation under Grant No. CCF-0634123. 

* Shanghai Key Laboratory of Trustworthy Computing, East China Normal University, 
Shanghai 200062, China, mwu@sei.ecnu.edu.cn, The author is supported in part by NSFC- 
44012140 and NSFC-90718041. 



1 



Systems 

3.1 Systems with Liouvillian Sequences as Solutions 

3.2 Normal Forms 

3.2.1 Normal Forms for General Systems 

3.2.2 Normal Forms for Systems of Prime Order . . . 

3.3 A Decision Procedure for Systems of Prime Order . . . 

3.3.1 The Decision Procedure for the Reducible Case . 

3.3.2 The Decision Procedure for the Irreducible Case 

3.4 Summary 



1 Introduction 

One of the initial and key applications of the Galois theory of linear dif- 
ferential equations is to characterize the solvability of such equations in 
terms of liouvillian functions, i.e., functions built up iteratively from ratio- 
nal functions using exponentiation, integration and algebraic functions ( |10l 
Appendix 6], |25l Chapters 1.5 and 4]). From the differential Galois theory, 
a linear differential equation 

L{y) = + a„_iy("-i) + • . • + aoy = 0, 

over C{x) can be solved in these terms if and only if its Galois group has a 
solvable subgroup of finite index. This allows us to conclude that if a linear 
differential equation L{y) = has a liouvillian solution then it has a solution 
of the form e^^ ^ where / is an algebraic function. This characterization is the 
foundation of many algorithms that allow one to decide if an equation has 
such solutions and find them if they exist. This theory and these algorithms 
have been developed for systems of matrix form Y' = AY as well as for more 
general coefficients. 

For the case of difference equations, the situation is in many ways not as 
well developed. A Galois theory of linear difference equations is developed 
in [26] . Later on in [13] , a notion of solving in liouvillian terms is introduced 
for linear difference equations of the form 

L{y) = y{x + n) + an-iy{x + n - I) -\ h aoy(a;) = with G C(x) 

and for difference equations in matrix form Y[x + 1) = AY{x) where A is a 
matrix over C{x). In [13], a characterization of solving in liouvillian terms 
is presented in terms of Galois groups and an algorithm is given to decide 



2 



whether a linear difference equation can be solved in liouvillian terms. Re- 
mark that in [13] , solutions are considered as equivalence classes of sequences 
of complex numbers y = (y(0), y{l), ■ ■ ■) where two sequences are equivalent 
if they agree from some point onward. A rational function / is identified 
with the sequence of the rational values (/(O), /(I), . . .). The addition and 
multiplication on sequences are defined elementwise. Liouvillian sequences 
are built up from rational sequences by successively adjoining solutions of 
the equations of the form y{x + 1) = a.{x)y{x) or y{x + 1) — y(x) = h{x) 
and using addition, multiplication and interlacing to define new sequences 
(the interlacing of u = {uq,ui, . . = {vq,vi,...) is (uq, fo, iti, fi, • • •))• 
Similar to the differential case, a linear difference equation L{y) = can 
be solved in terms of liouvillian sequences if and only if its Galois group 
has a solvable subgroup of finite index. When this is the case, L{y) = 
has a solution that is the interlacing of hypergeometric sequences, and [13] 
shows how to decide if this is the case. The paper |13j also gives examples 
of equations which have no hypergeometric solutions but do have solutions 
that are interlacings of hypergeometric solutions. Similar results also apply 
to difference equations in matrix form. 

We now consider systems of linear difference-differential equations of the form 
Y{x + l,t) =AY{x,t), ^^^^^ = BY{x,t) 

where A and B are square matrices over C(x,t) and A is invertible. In [TJ 
[21} [22| [28] , some theories and algorithms have been developed on determin- 
ing reducibility and existence of hyperexponential solutions of such systems. 
However, as in the pure difference case, there are systems which have no 
hyperexponential solutions but have solutions that are interlacings of hy- 
perexponential solutions. In this paper, we shall use a Galois theory that 
appears as a special case of the Galois theory developed in [12l Appendix] 
to characterize liouvillian solutions of linear difference- differential systems, 
and then devise an algorithm to determine whether such a system has liou- 
villian solutions when the order of the system is prime. 

Throughout the paper, we use (•)^ to denote the transpose of a vector or 
matrix and det(-) to denote the determinant of a square matrix. The sym- 
bols Z>o and Z>o represent the set of nonnegative integers and the set of 
positive integers, respectively. Denote by 1 the identity map on the sets in 
discussion. For a field k, denote by Qlnik) the set of all nxn matrices over k 
and by GL„(A:) the set of all n x n invertible matrices over k. All difference- 
differential systems of the form {a{Y) = AY, 6{Y) = BY} with A £ GL„(A;) 
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and B G Qln{k) that are in discussion in the paper are assumed to be inte- 
grable. 

The paper is organized as fohows. In Section [2l we wih first review some 
Galois theoretic results in [12] and then show that the Galois group of a 
linear difference-differential system is solvable by finite group if and only if 
a certain associated system has solutions in a tower built up using gener- 
alizations of liouvillian extensions. In Section [3l we show that irreducible 
systems with liouvillian solutions must be equivalent to systems of partic- 
ular form (Theorem 13. lip and refine this result for systems of prime order. 
We propose an algorithm for deciding if linear difference-differential systems 
of prime order have liouvillian solutions. At last, some examples are illus- 
trated. 

We would like to thank Reinhart Shaefke for supplying a simple proof of 
Lemma 13.251 

2 Galois Theory 

2.1 Picard-Vessiot extensions and Galois groups 

In [12], a general Galois theory is presented for linear integrable systems of 
difference-differential equations involving parameters. When there exists no 
parameters this theory yields a Galois theory of difference-differential sys- 
tems as above. Let us recall some notation and results in |12j . 

A ad-ring R \s a. commutative ring with unit endowed with an automor- 
phism a and a derivation 5 satisfying a6 = 5(t. R is called a a6 -field when R 
is a field. An element c of i? is called a constant if cr(c) = c and 5{c) = 0, 
i.e., it is a constant with respect both a and 6. The set of constants of R, 
denoted by R'^^, is a subring, and it is a subfield if i? is a field. 

In this section, unless specified otherwise, we always let fc be a aJ-field of 
characteristic zero and with an algebraically closed field of constants. 

Consider a system of the form 

a{Y) = AY, 6{Y) = BY (1) 

where A S GL„(A;), B G Qlnik) and y is a vector of unknowns of size n. 
The integer n is called the order of the system ([1]). A a6-nng R is called 
a a5-Picard-Vessiot extension, or a a6-P\ extension for short, of k for the 
system ([T]) if it satisfies the following conditions 
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(i) i? is a simple ad-nng; 

(ii) there exists Z £ GL.„(i?) such that cr{Z) = AZ and S{Z) = BZ; 

[in) R = k[Z, ^^1(^2 -j ], that is, R is generated by entries of Z and the inverse 
of the determinant of Z. 

Note that if the system ([1]) has a a6-P\ extension, the commutativity of a 
and 5 imphes 

a{B) = 5{A)A-^ +ABA-\ 

which is called the integrability conditions for the system ([1]). Conversely, if 
the system ([T]) satisfies the above integrability conditions and the constants 
of k are algebraically closed, it is shown in [7j and [121 Appendix] that cr6- 
PV extensions for ([T]) exist and are unique up to a5-k isomorphisms. 

The following notation will be used throughout the paper. 

Notation 2.1 Let A be a square matrix over a ad-ring. For a positive inte- 
ger m, denote Am = ct"^~ (A) ■ ■ ■ a{A)A. For a linear algebraic group G ^ G 
represents the identity component of G. 

Lemma 2.2 [Lemma 6.8 in [12]] Let k be a aS-field and R a simple a5-ring, 
finitely generated over k as a aS-ring. Then there are idempotents eg, ... , Cs-i 
in R such that 

(i) R = eoR® • • • e Cs-iR; 

(a) a permutes the set {cqR,--- ,es-iR}. Moreover, leaves each eiR 
invariant; 

[Hi) each CiR is a domain and a simple a^6-ring. 

The following lemma is an analogue to Lemma 1.26 in [25j . 

Lemma 2.3 Let k be ad- field, R be a ad-PV extension for the system 

a{Y) = AY, 5Y = BY with A G GL„(A;) and B G g/„(A:) 

and Co, ei, . . . , Cs-i be as in Lemma \2.2l Then each CiR is a a^6-PV exten- 
sion of k for the system {a^(Y) = AsY,6{Y) = BY}. 
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Proof. Let Rhea a6-PY extension for {a{Y)=AY,6{Y)=BY} and F be 
a fundamental matrix over R for the system. By Lemma 12.21 each Cji? is a 
simple o"*(5-ring. Clearly, eiF are the solutions of {a^{Y)=AsY,6(Y)=BY} 
since cr*(ei) = Cj and (^(cj) = for z = 0, . . . , s—1. Assume that Cj det(F) = 
for some i. Then 

a^{ei det(F)) = e^^. ^ det(^j) det(F) = 

and thus ej+j mod s det(F) = for j = 1, • • • , s — 1. Therefore 

det(F) = (eo + • • • + e,_i) det(F) = 0, 

a contradiction. So CiF is a fundamental matrix for {a'^{Y)=AsY, 5{Y)=BY} 
for each i. Moreover, CiR = k[eiF, ^-gj^-p^] for each i. This completes the 
proof. □ 

Corollary 2.4 Let d > 1 be a divisor of s. Suppose that there exist idem- 

potents Co, ... , eg-i in R such that R = e^R © eiR © • • • © Cg^iR. Then 

——1 

for i = 0, . . . ,s — 1, the subring (Bj^q Ci+jdR of R is a a'^6-PV extension 
of k for the system 

a\Y) = AdY, 6iY) = BY. 

Here we use a cyclic notation for the indices {0, • • • , s — 1}. 

Proof. The proof is similar to that of Lemma 12.31 □ 

Definition 2.5 Let A £ GL„(A;), B G Qln{k) andY be a vector of unknowns 
of size n. Let R be a a5-PV extension for {(j(Y) = AY, S{Y) = BY}. 
The group consisting of all a 5 -k- automorphisms of R is called the a5-Galois 
group for the system and denoted Gal{R/k). 

Denote by Gal{eoR/k) the o-^(5-Galois group for {a'{Y) = AY,6{Y) = BY}. 
Without loss of generality, we assume that cr(ej) = e^_^^ mod s- Construct a 
map r from Gal(eo-R/A:) to Gal(i?/A:) as follows. Let <p G Gal(eo-R/A;). For 
any r = tq + ri + • • • + r^-i G R with rj G ejR for j = 0, . . . , s — 1, define 

s-l 

rM(r) = ^a^ipa-^{rj). 

j=0 

Let (j) G Gal(i?/A;). Clearly, <j) permutes the e^'s by the proof of Lemma 6.8 
in [12]. Define a map A : Gal{R/k) Z/sZ to be A((/)) = i if (/>(eo) = e^. 
We then have the following 
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Lemma 2.6 Let R he a aS-PV extension for {(j{Y) = AY, 6{Y) = BY} 
where A G GLn{k) and B € Qlnik)- Let V and A he stated ahove. Then 
the map T is well-defined, i.e., if G Gal(eoii/A;) implies T{ip) G Gal(i?//c). 
Moreover, the sequence of group homomorphisms 

— > Gal{eoR/k) Gal{R/k) ^ Z/sZ — > 

is exact. 

Proof. The proof is similar to that of Corollary 1.17 in [25]. □ 

Lemma 2.7 Suppose that k has no proper algebraic a6-field extension and R 
is a a6-PV extension for {a{Y) = AY, 6{Y) = BY} where A G GL„(A;) 
and B e Qln{k). Then Gal{R/kf = Gal{eoR/k). 

Proof. Let k be the algebraic closure of k in the quotient field of cqR. 
Then Gal(eo-R/A;)° = Gal{eoR/k). Since k has no proper algebraic a5-G.eld 
extension, we have k = k and therefore Gal(eo-R/A:) = Gal(eo-R/A;)'^. Prom 
Lemma |2.6| it follows that Gal(eo-R/fe) is a closed subgroup of Gal{R/k) of 
finite index. The proposition in [n\ p. 53] then implies the lemma. □ 

Prom |12j . we know that a a5-FY extension R over k is the coordinate ring 
of a Gal(itI/A;)-torsor over k. Let E be an algebraically closed differential 
field with a derivation 5. Clearly, E{x) becomes a crj-field endowed with 
the extended derivation 5 such that 6{x) = and with an automorphism a 
on E{x) given by a\E = 1 and a{x) = x + 1. Por such a field E{x), we will 
get an analogue of Proposition 1.20 in [25]. Before we state the result, let 
us look at the following 

Lemma 2.8 Let k be a differential field with a derivation 5, S a differential 
ring extension of k and L a differential radical ideal of S. Suppose that S 
is Noetherian as an algebraic ring and that L has the minimal prime ideal 
decomposition n*^^Pj as an algebraic ideal. Then Pi are differential ideals 
fori = l,...,t. 

Proof. Let /i e Pi and select /j G Pj \ Pi for i = 2, • • • , t. Then / = 
/1/2 ■ ■ ■ ft ^ I- By taking a derivation on both sides, we have 

5(/) = 5(/i)/2---/t + M(/2)---/t + --- + /i/2---5(/t)e/, 

which implies (5(/i)/2 ■ • • /t G Pi- So G Pi and Pi is a differential ideal. 
The proofs for other Pj's are similar. □ 
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Now let S* be a finitely generated a^-ring over k and / a radical c^-ideal 
of S. Suppose that S is Noetherian as an algebraic ring and / = nf^j^Pj 
is the minimal prime decomposition of / as an algebraic ideal. Since S 
is Noetherian, we have = /, which implies that a permutes the Pj's. 
From Lemma [2 .St each Pj is a differential ideal. Therefore if {Pi}i^j with J 
a subset of {1, ... ,s} is left invariant under the action of u, then Di^jPi 
is a (T(5-ideal. We then have the following result. We will use the following 
notation: if V is a variety defiend over a ring ko and ki is a ring containing /cq, 
we denote by V{ki) the points of V with coordinates in ki. 

Proposition 2.9 Let k = E{x) be as in the paragraph preceding Lemma \2.8l 
R a a5-PV extension of k for the system {o'{Y) = AY, 6{Y) = BY} 
where A £ GL„(A;) and B G Qln{k), and G = Gal[R/k). Then the corre- 
sponding G-torsor Z has a point which is rational over k and Z[k) and G{k) 
are isomorphic. Moreover, G/G^ is cyclic. 

Proof. The notation and proof will follow that of Proposition 1.20 in [25j . 
Let Zq, ■ ■ ■ ,Zt-i be the ^-components of Z. By Lemma |2.8| the defining 
ideals Pi of Zi are differential ideals. As in the proof of Proposition 1.20 
in [25], there exists B £ ZQ{k) such that Zq = BG^ and Z = BGj^ where G^ 

denotes the variety G over k. Since Z[k) is r-invariant, we have 

BG~^ = t{BGj^) = A-^a{B)G-„ 

which implies ^-^-^(P) G G{k). There exists N G Gik"^) such that 

B-^A-^a{B) G G°{k)N. 

Let H be the group generated by G^ and A^. One sees that t{BHj^ = BHj^ 
and therefore the defining ideal / of BHj, is cr- invariant. Since the set BH^ 
is the union of some of the Zj, / is of the form rijgjPj with J a subset 
of {0, 1, • • • , t— 1}. Hence / is a crS-ideal because each Pi is a differential ideal. 
Since the defining ideal / of Z is a maximal ad-ideal, it follows that 1 = 1 
and so H = G. □ 

From B~^A~^a{B) G G^{k)N in the proof of Proposition 12.91 we conclude 
that A is a generator of the cyclic group G/G^ . 

Definition 2.10 Let k he a a5-field and R he a ct5-PV extension of k for 
the system 

a{Y) = AY, 5{Y) = BY 
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with A G GLn{k) and B G g[„(A;), and V C i?" he the solution space of the 
system. The system is said to be irreducible over k if V has no nontriv- 
ial Gal{R/k) -invariant subspaces. 



In a manner similar to the purely differential case [261 p. 56], one can show 
that a system {a{Y) = AY, d{Y) = BY} with A G GL„(A;) and B G 0[„(fc) 
is reducible over k if an only if there exists U G GL„(/c) such that a change 
of variables Z = UY yields a{Z) = AZ, 5{Z) = BZ with 



We end this section with a concrete way of realizing a cJ-Picard-Vessiot 
extensions for linear difference-differential systems over fields of particular 
form. We proceed in a manner similar to that of [261 Chapter 3]. 

Let be a differential field with a derivation b. Denote by Sk the set of all 
sequences of the form a = (oq, oi, . . .) with a(i) = ai £ K. Define an equiv- 
alence relation on Sk as follows: any two sequences a and b are equivalent 
if there exists N G Z>o such that a(n) = b(n) for all n > N. Denote by Sk 
the set of equivalence classes of Sk modulo the equivalence relation. One 
sees that Sk forms a differential ring with the addition, multiplication and 
a derivation 6 defined on Sk coordinatewise. Clearly, the map a given by 
(7((ao, oi, . . .)) = (ai, a2, . . .) is an automorphism Sk that commutes with 
the derivation 6. In addition, any element e € K is identified with (e, e, . . .). 
So we can regard as a (differential) subfield of Sk ■ 

Let E be an algebraically closed differential field with a derivation 5. Con- 
struct an automorphism a on E{x) given by a\E = 1 and 0"(x) = x + 1 
and extend 5 to be a derivation 6 on E{x) such that 6{x) = 0. Assume 
that K is a differential field extension of E with an extended derivation 6. 
The map E{x) Sk given by / i-^ (0, . . . , 0, /(TV), /(TV + 1), . . .), where N 
is a non-negative integer such that / has no poles at integers > A^, induces 
a (T(5-embedding of E{x) into Sk- Consequently, we may identify any ma- 
trix M over E{x) with a sequence of matrices (0, . . . , 0, M(A^), M(A^+1), • • • ) 
where M{i) means the evaluation of the entries of M at x = i. So we have 
the following 

Proposition 2.11 Let E C K, E{x), Sk be as above and let k = E{x). 
Assume that E and K have the same algebraically closed field of constants 




Note that A and B are again n x n matrices over k. 



9 



as differential fields. Let R be a a5-PV extension for the system 

a{Y) = AY, 5{Y) = BY 

where A G GL„(A;) and B G glnik). Let N G Z>o such that A{m) and B(m) 
are defined and det(A(m)) 7^ for all m > N and assume that S{Y) = 
B{N)Y has a fundamental matrix Z £ GL„(i^). Then there exists a a5- 
k-monomorphism of R into Sk- Moreover, the entries of any solution 
of {a{Y) = AY, 6{Y) = BY} in lies in the image of R in Sk- 

Proof. Let R = k\Y, dcty ]/-^ be the cj5-PV ring extension for the sys- 
tem {a{Y) = AY,5{Y) = BY} and let G be its Galois group. From 
Proposition 12.91 the corresponding torsor has a point P with coordinates 
in k. This implies that if we introduce a new matrix of variables X and 
let Y = PX, then R = k[X, -^^j^]/J where J is the defining ideal of G. 
Furthermore, a{X) = AX and 5{X) = BX where A = a{P)-^AP G G{k) 
and B = P~^BP — P~^5{P) G Q{k) where Q is the lie algebra of G. 

Define recursively a sequence of matrices G GL.„(/^) for m > N: 

Zfq = Z and Zm+i = A{m)Zm for any m > N. 

The integrability condition on A and B implies that Zm satisfies S{Y) = 
B{m)Y for any m > N and so Z = (. . . , Z^, -^Tv+i, • • •) is a fundamental 
matrix in GL„(cSa') of {a{Y) = AY, 5{Y) = BY). 

Remark that the key to proving the proposition is to show that Z generates 
a (jJ-Picard-Vessiot extension. Unfortunately, we do not see a direct way to 
show this and our proof is a little circuitous. Clearly, U := P~^Z satifies 
that a{U) = AU and d{U) = BU. Then 6{U{N')) = B{N')U{N') for a 
sufficiently large A^' and therefore K contains a (differential) Picard-Vessiot 
extension of E for the equation d{U) = BU. Since B{N') G Q{k), Propo- 
sition 1.31 (and its proof) in [26] together with the uniqueness of Picard- 
Vessiot extensions imply that there exists V G G{K) such that 5{y) = 
B{N')V. Define V e Sr hy V{N') = V and V{m + 1) = A{m)V{m) 
for m > N'. Then a{V) = AV and 6{V) = BV and V G G{Sk). This 
implies that the map from R = k[X, -^^j^]/ J to Sk given by X 1— > F is 
a (T(5-A;-homomorphism. Since / is a maximal o"5-ideal, this map must be 
injective, and so is the desired embedding from R into Sk- 

Let Ty G 5^ be a solution of {a{Y) = AY, 5{Y) = BY}. For a suffi- 
ciently large M, W{M) is defined and is a solution of 5{Y) = B{M)Y 
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and W{m + 1) = A{m)W{m) for m> M. Therefore W{M) = V{M)D for 
some constant vector D and thus W = VD G i?". It fohows that Z = PU 
also generates a crJ-Picard-Vessiot extension, as clahned in the above. □ 

Remarks 2.12 If K is a maximal Picard-Vessiot extension of E with the 
same constants (Zorn's lemma guarantees that such fields exist), then the 
hypothesis on the existence of Z in Proposition \2.11\ is always satisfied. 
Therefore for such a field K , Sk contains a a 5 -Picard-Vessiot ring for any 
system {a{Y) = AY, 5{Y) = BY}. 

2.2 Liouvillian solutions 

The Galois theory for linear differential equations is stated in terms of dif- 
ferential integral domains and fields [261 Chapter 1] and both theory and 
algorithms for finding liouvillian solutions are well developed |261 Chapters 
1.5, 4.1 - 4.4]. The main result is that the associated Picard-Vessiot exten- 
sion lies in a tower of fields built up by successively adjoining, exponentials, 
integrals and algebraics if and only if the associated Galois group has a solv- 
able identity component. For linear difference equations, the Galois theory 
is stated in terms of reduced rings and total rings of fractions. A general 
theory of liouvillian solutions has not been developed in the difference case. 
However, a case has been investigated in [13] where the coefficient field is of 
the form C{x) with a shift operator a : x ^ x + 1 and a\c = 1. In this sit- 
uation, solutions of linear difference equations are identified with sequences 
whose entries are in C. One says that a linear difference equation is solvable 
in terms of liouvillian sequences if it has a full set of solutions in a ring 
of sequences built up by successively adjoining to C sequences representing 
indefinite sums, indefinite products and inter lacings of previously defined 
sequences. The main result is that a linear difference equation is solvable in 
terms of liouvillian sequences if and only if its Galois group has a solvable 
identity component. 

In this paper we will combine the approaches for differential and difference 
cases to investigate the solvability of systems of mixed linear difference- 
differential equations over E{x) where E will always be an algebraically 
closed differential field unless specified otherwise, a{x) = x + 1 and a\E = 1. 

In this section, we will give a characterization of Galois groups for mixed 
difference-differential systems to have solvable identity component in terms 
of liouvillian towers over an arbitrary fT(5-field k with algebraically closed 
constants. Then we will define a notion of liouvillian sequences and show 
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that having a fuU set of solutions of this type imphes that the Galois group 
has solvable identity component (Proposition 12.23]) . In a later result (Propo- 
sition [STT]), we will show the converse is true as well. 



Liouvillian extensions for aS -fields are defined in the usual way. 

Definition 2.13 Let k be a a5-field. A a6-field extension K of k is said to 
be liouvillian if there is a chain of a6- field extensions 



such that k = K" , i.e., K shares the same set of constants with k, 
and -fCj+i = Ki[ti) for i = 0, . . . , m — 1 where 

(1) tj is algebraic over Ki, or 

(2) a{ti) = riti and 5{ti) = r2ti withr\,r2 G Ki i.e., ti is hyperexponential 
over Ki, or 

(3) a{ti) -ti£ Ki and 6{ti) G Ki. 

We now define liouvillian solutions of mixed difference-differential systems. 
In the sequel, let k he a uJ-field with algebraically closed constants, R 
be a a5- PV extension of k for the system {(t(Y) = AY, d(Y) = BY} 
with A G GL.„(/c) and B G Qln{k)- Suppose that R has a decomposition 



and J^Q is the quotient field of cqR. Then J^q is a a^6-&eld. 

Definition 2.14 Let v = T^tZoVi be a solution of {a{Y)=AY, 5{Y)=BY} 
in i?" where Vi := CiV G CiR^. We say that v is liouvillian if the entries 
of vq lie in a 5 -liouvillian extension of k containing Tq. We say that the 
original system is solvable in liouvillian terms if each solution is liouvillian. 

Suppose that v is a solution of {(jiY) = AY, 5{Y) = BY} in i?". From 



A; = i^o C i^i C • • • C Km= K 



R = eoR ® eiR © • • • Cs-iR 



cr{v) = Av and 5{v) = Bv, 



it follows that 



a{vs-i) © o-(vo) © • • • © (t{Vs-2) 
S{vo) ® 6{vi) ® ■ ■ ■ ® 6{vs-i) 



Avo © Avi © • • • © Avs-i 
Bvo © Bvi © • • • © Bvs-i 
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which imphes that a{vi) = ^f^j+imods ^{"^i) — for i = 0, . . . , s—1. 
Hence cr^{v()) = AgVo and 6{vo) = Bvq, i.e., vq is a solution of the sys- 
tem {(T*(y) = AgY, d{Y) = BY}. Conversely, assume that vq is a solu- 
tion of {a'{Y) = AsY, 6{Y) = BY} in eoi?". Let Vi = A-^a{vi-i) for 
i = 1, ■ ■ ■ , s — 1 and v = vq(B - ■ ■ ® Wg-i- We then have cr{v) = Av. From the 
fact that cr(B) — 6{A)A~^ = ABA~^, one can easily check that 5{v) = Bv. 
Hence f is a solution of {(j{Y) = AY, d{Y) = BY} in i?". Moreover we have 
the following 

Proposition 2.15 The system {(7{Y) = AY, 6{Y) = BY} is solvable in 
liouvillian terms if and only if the system {a'^(Y) = AgY, d{Y) = BY} is 
solvable in liouvillian terms. 

Proof. Let Vi,--- ,Vn be a basis of the solution space for the system 
{(y{Y) = AY, 5{Y) = BY}. Since Vi is liouvillian, Vio = e^Vi is liouvillian 
for each i. It then suffices to show that ^lo, • • • , Vno are linearly independent 
over h"^ . Assume that there exist ci, • • • ,Cn e k''^ not ah zero, such that 

ciVu) H h c„y„o = 0. Letting Vn = eiVi we have Vn = A~^a{Vio) for i = 

1, . . . , n. Remark that k^^ = k^"^ since k"^ is algebraically closed. Therefore 

ciFii + • • • + CnVni = A-^a{ciViQ + ■■■ + c„T/„o) = 0. 

Similarly, cieiVi + • • • + CneiVn = for each i. Hence ciVi + • • • + CnVn = 0, 
a contradiction. 

Conversely, suppose that Vio, • • • , V^o is a basis of the solution space for the 
system {a^{Y) = AgY, 5{Y) = BY} and that all the Vi^s are liouvillian. 
For i = 1, • • • ,n and k = 1, - ■ ■ , s — 1, let 

Vik = A~'^a{Vik-i) and V, = Vio®Vii® ■ ■ ■ ®Vi^s-i- 

Then each Vi is a solution of {cj{Y) = AY, 5{Y) = BY}. Clearly, Vi, • • • , K 
are linearly independent over k"^ . This concludes the proposition. □ 

Theorem 2.16 The system {(t(Y) = AY, 6{Y) = BY} is solvable in liou- 
villian terms if and only if Gal{R/k)^ is solvable. 

Proof. By Proposition 12.151 {a{Y) = AY, 6{Y) = BY} is solvable in liou- 
villian terms if and only if the associated system {a'^{Y) = AgY, 5{Y)=BY} 
is solvable in liouvillian terms. By Lemma 12.61 Gal{eQR/k) is a subgroup 
of Gal{R/k) of finite index. Then Gal(R/kf = G&l{eoR/kf. Hence it suf- 
fices to show that the associated system is solvable in liouvillian terms if and 
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only if Gal(eoi?/A;)'^ is solvable. Note that the PV extension of k for the as- 
sociated system is a domain. Thus the proof is similar to that in differential 
case. □ 

Let E he a differential field with algebraically closed constants and L be a 
field containing E which satisfies the following conditions: 

(i) L is a differential field extension of E having the same field of constants 
as E; 

(ii) every element in L is liouvillian over E; 

(iii) L is maximal with respect to (i) and (ii). 

Zorn's Lemma guarantees that such a field exists. We refer to L as a maximal 
liouvillian extension of E. One can show that any differential liouvillian 
extension of E having the same field of constants as E can be embedded 
into a maximal liouvillian extension and that any two maximal liouvillian 
extensions of E are isomorphic over E as differential fields. 

Recall that for two sequences a and b, and for a nonnegative integer m, b 
is called the ith m-inter lacing ( [13^ Definition 3.2]) of a with zeroes if 

h{mn + i) = a(n) and b(r) = for any r ^ i mod m. 

A sequence a is called the ith m-section of b if a(mn + i) = h{mn + i) and 
a(r) = for r ^ i mod m. 

We now turn to a definition of liouvillian sequences. 

Definition 2.17 Let E be a differential field with algebraically closed con- 
stants, a be an automorphism on E{x) satisfying cr{x) = x + 1 and a\E = 1- 
Let L be a maximal liouvillian extension of E. The ring of liouvillian se- 
quences over E{x) is the smallest subring C of Sl such that 

(1) L{x) C C; 

(2) For a.e Sl, a G £ if and only if (j(a) G C; 

(3) Supposing cr{h) = ah with a, b G Sl, then a G E{x) implies b G £. b 
is called a hypergeometric sequence over E[x); 

(4) Supposing cr{h) = a -|- b with a, b G 5^, then a G £ implies b G £; 

(5) For a G Sl, a G £ implies that b G £, where b is the ith in-interlacing 
of a with zeroes for some m G Z>o and < i < m — 1. 
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Set k = E{x) and let C be the ring of liouvillian sequences over k. From the 
remarks in [131 p. 243], if b € 5/, belongs to C then the ith m-section of b 
also belongs to C for any i and m. We claim that £ is a ad-rvag. Since C 
can be constructed inductively using (1) - (5) above, it is enough to show 
the following statements: 

(1') If / G L{x) then 5{f) E L(x); 

(2') If a G Sl, a, (5(a) G C then 5(a(a)) G £; 

(3') Supposing (j(b) = ab with a, b G 5l, then a, 5(a) G £'(x) implies 
8{h) G £; 

(4') Supposing o"(b) = a + b with a, b G Sl, then a, (5(a) G C implies 
<5(b) G £; 

(5') For a G Sl-, a, (5(a) G C implies that (5(b) G £, where b is the ith 
m-inter lacing of a with zeroes for some m G Z>o and < i < m — 1. 

Verifying (!'), (2') an (5') is straighfoward. To verify (3'), suppose that b G £ 
with (T(b) = ab. Set y = (5(b). Then cr(y) — ay = (5(a)b. Since a G fe, b is 
invertible. Then u := (5(b)/b satisfies cr{u) — u = (5(a)/a. Therefore u G £ 
and so 6(b) = ub G C. To verify (4'), suppose that b G £ with <7(b) = a+b. 
Then cr{6{h)) = (5(b) + (5(a) which implies that (5(b) G £. 

A vector is said to be hypergeonietric over k if it can be written as Wh 
where W £ k^ and /i is a hypergeometric sequence over k. For any positive 
integer d, we can construct a solution of {a'^{Y) = A^Y, 6{Y) = BY^ by 
interlacing as indicated below. Consider a set of new systems 

S,-: {a{Z) = Ad{dx+i)Z, 5{Z) = B{dx + j)Z}, j = 0,---,d-l, 

where A(i{dx + j) and B{dx + j) mean replacing x hy dx + j in each entry 
of Ad and B, respectively. Clearly, these new systems are all integrable. 
Moreover, we have the following 

Proposition 2.18 Let Vj be a solution of Sj in where Sk is as in 
Remark \2.1S\ and Wj be the jth d-interlacing with zeros of Vj for j = 
0, • • • , d — 1. Then 

W = Wo + --- + Wd-i 
is a solution of {a'^{Y) = AdY, 5{Y) = BY]. 
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Proof. It suffices to show that is a solution of {cj^(y)=^dy, 6{Y)=BY} 
for each j. Let Vj = {Vj{0),Vj{l), ■ ■ ■ ) and Wj = {Wj{0),Wj{l), ■ ■ ■ ) for 
j = 0,--- ,d-l. Then 

V,{i + 1) = Ad{di+j)V,{i), S{Vj{^) = B{di+j)Vj{i) 

and therefore by definition of interlacing, for i > and I = - ■ ■ ,d— 1, 

Wj{di +j + d) = Vj{i + 1) = Ad{di + j)Vj{i) = Ad{di + j)Wj{di + i), 
+j + (i + £) = = + + j + 1). 

So cr'^(W;,) = AdVFj. It is clear that 5{Wj) = BWj. So the proposition 



Definition 2.19 Let E he an algebraically closed differential field with a 
derivation 5. Suppose that 5 is extended to be a derivation 6 on E{x) 
such that 5{x) = and a is an automorphism on E{x) such that a\E = 1 
and a{x) = x + 1. A system 



over E{x) is said to be solvable in terms of liouvillian sequences if the a5-PV 
extension of this system embeds, over E{x), into C, the ring of liouvillian 
sequences over E{x). 

Clearly, Definition 12.191 generalizes the notion of solvability in liouvillian 
terms for linear differential equations and that for linear difference equations. 
In addition, we shall show in Proposition 12.231 that this property is also 
equivalent to the identity component of the Galois group being solvable. 
However this property is not equivalent to having a solution in the ring Sl, 
as shown by the following 

Example 2.20 The Hermite polynomials 



satisfy a linear differential equation with respect to t and a difference equa- 
tion with respect to n (cf., JBil, Ch. 10.13). In matrix terms, the vector 
Y{n,t) = {H{n,t),H{n + l,t))'^ satisfies 



holds. □ 



a{Y) = AY, 5{Y) = BY 




[n/2] 
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so this system has a solution in Sl and one can use the techniques of Sec- 
tion \3.3\ to show that this system have no solutions in C. 

Concerning total rings of fractions, we have the following 

Lemma 2.21 Let k be a ad-field and S a a6-PV extension ofk. A nonzero 
element r £ S is a zero divisor in S if and only if there exists j G Z>o such 

that nio^'(o = o. 

Proof. Suppose that Ylj-Qcr'^ir) = for some j G Z>o and let j be min- 
imal with respect to this assumption. If = rii=i'^*('^) ~ '^(TiiZo ^ 
then we would also have nto'^H'') — contradicting the minimality of 
j. Therefore r is a zero divisor. Now suppose that r is a zero divisor. 
We write S = (B{ZoSi where the Si are domains and o-{Si) = -Sj+imodj- 
Then r = "^{Zq ri with G Si. Writing r = (ro, . . . , J^j-i), one sees that r 
is a zero divisor in S if and only if some is zero. Assume ro = 0. Then 

cT(r) = (a(rj„i),0,...), a2(r) = (a2(rj_2),fT'(rj-i),0,...), ... 

soO = n£o^'(0- □ 
Consequently, we have 

Corollary 2.22 Let k and S be as above and let S be a a-subring of S. An 
element r £ S is a zero divisor in S if and only if it is a zero divisor in S. 
Therefore the total ring of fractions of S embeds in the total ring of fractions 
ofS. 

Proposition 2.23 Let E[x) be as in Definition \2.19[ If a system 

a{Y) = AY, 6{Y) = BY 

over E[x) is solvable in terms of liouvillian sequences, then the identity 
component of its Galois group is solvable. 

Proof. Let R, L and C be the (T(5-PV ring, a maximally liouvillian extension 
of E and the ring of liouvillian sequences for the given system, respectively. 
Then R C C Consider the following diagram 

RL{x) 

I \ 
R L{x) 
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\ / 

RnL{x) 



E{x) 

where RL{x) is the field generated by R and L(x). We will show that 

(i) RL{x) C £ is a a6-PY extension of L{x) and its Galois group has a 
solvable identity component. 

(ii) The crJ-Galois group of RL{x) over L{x) is isomorphic to the sub- 
group H of the Galois group G oi R over E[x) that leaves the quotient 
field of Rr\L{x) fixed. Moreover, H \s a. closed normal subgroup of G. 

(iii) G/H is the Galois group of the quotient field of n L{x) over E{x) 
and has solvable identity component. 

Once the above claims are proven, the group G has a solvable identity com- 
ponent since both H and G/H have solvable identity component. 

To prove (i), we consider {(j(Y) = AY, 5{Y) = BY} as a system over L{x). 
Since R C C, then for a sufficiently large N there is a fundamental ma- 
trix of S{Y) = B{N)Y with entries in L. Applying Proposition 12. IH 
we conclude that Sl contains a crJ-PV extension T of L{x) for the sys- 
tem {a{Y) = AY, 5{Y) = BY} and that RcT. This imphes that RL{x) 
is a o"(5-PV extension of L{x). Proposition 4.1 of [25] implies that RL{x) 
is also a difference Picard-Vessiot extension for criY) = AY and the results 
of [H] then imply that the difference Galois group has solvable identity 
component. The (T(5-Galois group is a subgroup of this latter group and its 
identity component is a subgroup of the identity component of the larger 
group. Therefore it is also solvable. 

To prove (ii), let Fl be the total ring of fractions of RL{x) and Fe the 
total ring of fractions of R. Corollary 12.221 implies that we can regard Fe 
as a subset of Fl. The elements of the cr 5-Galois group G of F^ over L{x) 
restrict to automorphisms of Fe over E{x), and this gives a homomorphism 
of this group into H. Clearly the image is closed and the set of elements left 
fixed by this group is L{x) n Fe, the quotient field of L[x) H R. Therefore 
this image is H. 

Before proving that H is normal in G, we first show that 

FEDLix) = {FEnL){x). 
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Since x G Fe we have FEr\L{x) D (F£;nL)(x). To get the reverse inclusion, 
let / G Fe n L{x) and write 



OrX"" + . . . + ao , ^ r 
bsX^ + ... + bo 

where the numerator and denominator are relatively prime. We then have 
that {x*/, ...,/, x'', 1} are linearly dependent over L. Since L 

is the set of u-invariant elements of L(x), the cassoratian of these elements 
must vanish ([8], p. 271). This implies further that these elements are linearly 
dependent over the cj-invariant elements of the field Fe D L[x). Therefore 
there exist cr- invariant elements cij, hj G Fe H L for i = 0, 1, . . . , r and j = 
0, 1, . . . , s, not all zero, such that 

arx"" H h ao - {bsx'' f H h hof) = 0. 

Since x is transcendental over cj-invariant elements, there exists at least one 
hi which is not zero. Hence 

bgX^ + ... + bo 



To show that H is normal, it now suffices to prove that any cr5-automorphism 
of Fe over k leaves the field FEr\L{x) = {FEr\L){x) invariant. Note that L 
is the set of cr-invariant elements of F^ and so Fe H L is the set of a-invariant 
elements of Fe- This set is clearly preserved by any c7(5-automorphism. 

To prove (iii), note that since H is normal, the field {Fe H L){x) is a a6-P\ 
extension of E{x). Furthermore, {Fe H L){x) lies in the liouvillian exten- 
sion L{x) of E{x). Theorem 12.161 implies that its Galois group is solvable. □ 

Later in Proposition 13.71 we will show that the converse is true as well. 
2.3 From e^R to R 

Throughout this section, let E be an algebraically closed differential field 
with a derivation 5 whose extension on E{x) satisfies 6{x) = 0, and let a be 
an automorphism on E{x) such that a\E = 1 and (j{x) = x+\. Set k = E{x) 
and let C be the ring of liouvillian sequences over k. In this section, we shall 
show how to construct a solution in £" of {(t(Y) = AY., ^iX) = BY} from a 
solution in £" of its associated system {a'^iY) = AciY,5{Y) = BY}. More- 
over, we shall prove that if the associated system is equivalent over /c to a 
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diagonal form, so is the original system under the assumption that A is of 
particular form. 

Let F G be a nonzero solution of {o-^(y) = AdY,6(Y) = BY} and N £ 
Z>o be such that V{N) ^ 0, A{j) and B{j) are weh defined and det(A(j)) / 
for j > N. We define a vector W in the following way: 

W{N) = ViN) and W{j + 1) = A{j)W{j) for j > N. 

Since 6{V{N)) = B{N)V{N), the integrability condition on a and 6 implies 
that is a nonzero solution of {cr{Y) = AY, S{Y) = BY}. The proposition 
below says that W is also in 

Proposition 2.24 Let W be as above. Then W G C"-. 

Proof. Let N = M + m where ^, m G Z>o, < m < £ and Vq be the mth d- 
section of V. Then Vb is a solution of {o-'^(y) = AdY,5{Y) = BY} and 
Vq£ C Let 

^i(j) = + 1) for J > iV and i = 1, • • • , d - 1, 

<mdU = Vq + Vi + ■■■ + Vd-i. Then U £ We shah prove that W = U. 
Note that for j > iV, 

^i(j) = A{j)~'A{j + . . . A(j + i - 1)- Vo(j + i). 

In particular, Vd_i(j + 1) = ^(j)Vo(j). Then ^^(iV) = for i = 1, • • • , d- 1. 
Therefore VF(7V) = l/o(A^) + Vi{N) + --- + Vd-i{N) = U{N) and 

U{j + 1) = Vo(j + 1) + + !) + ••• + Vd-i{j + 1) 

= AimU) + ^(j)y2(j) + • • • + ^(j)^o(j) = A{j)U{j) 

for j > N. Hence TV = C/ G □ 

Lemma 2.25 Assume that w £ k satisfies 

a'{w) = a'~\b) ■ ■ ■ a{b)bw 

where s £ Z>o, b £ k\ {0} and b = x'^ + bix"^"^ + • • • with v £'L and hi £ E. 
Then a{w) = bw. 
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Proof. We have 

(^) = (^^"'(^) • • • = <^'^Hb) ■ ■ ■ a{b)b^. 

Since G k, = cw for some c € k'^^ = k^ . Note that w and h are 
rational functions in x. Expanding w and b as Laurent series at x = oo. By 
comparing the coefficients, we get c = 1, so a{w) = bw. □ 

Theorem 2.26 Let A = diag(ai, • • • , a„) where Oj G ^\{0} and = cx^^ + 
aiix^^~^ + • • • tfii/i G Z and c, G E. Assume that the system {a'^iY) = 
AfiY,6{Y) = BY} is equivalent over k to 

a''(Y) = AaY, <5(y) = diag(5i, • • • , 

where bi £ k for i = 1, . . . , n. Then {(t{Y) = AY, d{Y) = BY} is equivalent 
over k to 

a{Y) = AY, 6(Y) = diag(6i, • • • , bn)Y. 
Proof. From the assumption, there exists G G GL„(/c) such that 
a''{G)Ad = AdG, G-'BG - G~'S{G) = diag(5i, • • • , b^). 
It then suffices to prove that a{G)A = AG. Let G = {gij)nxn- Then 

{(^'^{gu) - 9ii = 0, i = l,--- ,n; 

a%,^) = a-^-i (l) . . . a (^) ^g,„ 1 < z ^ j < n. 

By Lemma I2.25| cr{gij) = ^gij for all i,j = 1, • • ■ ,n. This imphes that 
a{G)A = AG. a 

2.4 From C{x, t) to C(t)(a;) 

In the following sections, we always let fco be the ad-Held C(t, x) with an 
automorphism a : x i— > x+1 and a derivation <5 = and let k be its extension 
field C(t)(x). Consider a system of difference-differential equations 

a{Y) = AY, 6{Y) = BY 

over ko where A G GL„(A:o) and B G Qlniko)- We shall analyze this sys- 
tem by focusing on its difference part ^{Y) = AY and use techniques from 
the theory of difference equations. In this latter theory, one assumes that 



21 



the fixed field of a, that is, the cj-constants, are algebraically closed. For 
this reason we will need to consider properties of (j(Y) = AY over k as 
well as over /cq- We shall first show that the cJ-Galois group of this sys- 
tem over k can be identified with a normal subgroup of the Galois group 
of the same system over kQ, and then conclude some results on orders of 
the factors. For example, if the above system is irreducible over /cq, it is 
possible that the system is reducible over k. In this case, we will prove that 
the factors of the above system over k have the same order. A similar re- 
sult is well known for differential equations-if one makes a normal algebraic 
extension of the base field then the differential Galois group over this new 
field is a normal subgroup of the differential Galois group over the original 
field and an irreducible equation factors into factors of equal order. In the 
mixed difference-differential case or even the difference case, the fact that 
Picard-Vessiot extensions may contain zero divisors introduces some small 
complication. 

We start with some lemmas. Let R and Rq be the a5-FY extensions of k 
and kQ for the system {cr{Y) = AY,6{Y) = BY} respectively. 

Lemma 2.27 (i) There is a ko-monomorphism of ad-rings from Rq to R. 

(ii) Identify Rq with a suhring of R as in the first assertion. Suppose 
that Rq = foRo © • • • © fd^iRo where /ii?o is a domain, ff = fi 
and (j{fi) = fi+imodd, o-nd that R = e^R © • • • © Cg-iR is a similar 
decomposition of R. Then s = md for some m G Z>o. Moreover, after 
a possible renumbering of the fi, we have 

fi = ei + Ci+d H \- ei+(m-i)d for i = {),■■■ ,d-l. 

Proof, (i) Clearly, the ring Rq k becomes a a6-nng endowed with the 
actions (5(r©/i) = 5(r)©/i+r©(5(/i) and a{r®h) = a{r)®a{h) for any r £ Rq 
and h £ k. Since ko is a field, the two canonical embeddings 

Rq Rq ©fco 1 C i?o ©fco k and A; ^ 1 ©^^ k C Rq ©fc^ k 

are both injective, and clearly are homomorphisms of aJ-rings. Let M 
be a maximal o"(5-ideal of Rq ©^^ k and consider the ring (Rq ©^^ k)/M. 
Because Rq and k are both simple uJ-rings, the above embeddings factor 
through to (-Ro©fco k)/M and are still injective. Note that (-Ro©fco k)/M is a 
(T(5-PV extension of k for {(t{Y) = AY,5{Y) = BY}. So by uniqueness, we 
may write R = {Rq ©^q k)/M. Assume that Rq = kQ[Z, ^ ] where Z is a 
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fundamental matrix of the system. Let Z = Z mod M. One sees that Z is 
still a fundamental matrix of the system and that det(Z) 7^ 0. Hence 



R = ih 



z. ' 



' det(Z) 



(^ko k] /M = k 



z ' 



' det(Z) 



is a <t5-PV ring for the given system over k and Rq can be embedded into R. 

(ii) Write /o = X^^IoOjej with aj G Cji?. Squaring both sides yields /o = 
^^~Qa|ej, thus ajej = ajej. Since ejR is a domain, is either ej or 
for each j. The same holds for other /j's. Then for any i = 0, - ■ ■ ,d— 1, 
there is a subset Tj C {0, ••• ,s — 1} such that fi = '^Zji^T^^i- Assume 
that n Tj^ is not empty for two different and ii. Let / G Tj^ n Tj^. 
Since Eto /i = E^^o e,' = 1, 

d-l s-1 d-l s-1 

i=o j=o 1=0 jeTi j=o 

where p > and is the sum of all the e^'s with q ^ I. Multiplying 
both sides of the above equality by e/, we get pei = 0, a contradiction. 
Hence the Tj's form a partition of {0, • • • ,s— 1}. Since o"(/j) = /i+imodd 
and (j{ej) = ej+imods; one sees that the sets Tj have the same size and that 
a renumbering yields the conclusion. □ 

According to Lemma 12.271 we will always consider Rq as a subring of R 
and assume R = k[Z, -^^^^j in the sequel. In particular, we can view 

Lemma 2.28 Let^: Gsl^R/k) ^ Gal(-Ro/^o) be a map given by j{(j))=4>\Ro 
for any (j) € Gal(ii/A;). Then ^ is a monomorphism. Moreover, we can view 
the identity component of Gal{R/k) as a subgroup of that 0/ Gal(i?o//co)- 

Proof. Assume that R = k[Z, 5^]. Let (p £ Gal{R/k). If (t){Z) = Z[(t)]z 
for some [c/)]^ G GL„(C), then ^{4>){Z) = Z[(j)]2- Hence 7((/>) is an automor- 
phism of Rq over ko, that is, 7(0) G Gal(i?o//Eo)- Note that det(Z) 7^ 0. 
If 7(1^) = 1, then [(j)]2 = In, which implies that 0=1. So 7 is an in- 
jective homomorphism. Therefore, we can view Gsd{R/k) as a subgroup 
of Gal(i?o/A;o)- Since 7 is continuous in the Zariski topology, 7(Gal(i?/A;)'^) 
is in Gal(i?o/A:o)'^- So the lemma holds. □ 

Lemma 2.29 Gal{R/k)° = Gal(i2o/A:o)°- 
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Proof. Let G = Gal{R/k) and Go = Gal(i?o/A;o). From [l2], the a6-PY 
ring R (resp. Rq) is the coordinate ring of a G-torsor (resp. Go-torsor). 
From [201 p. 40, (3)], the Krull dimension of R (resp. Rq) equals the Krull 
dimension of G (resp. Go). Since all the components of a linear algebraic 
group are isomorphic as varieties, one sees that the Krull dimension of G 
(resp. Go) equals the Krull dimension of (resp. Gg). Since R is generated 
over Rq by the elements of k, by Proposition 2.2 and Corollary 2.3 in [20L 
p. 44], R is an integral ring extension of Rq. By [201 Corollary 2.13], the 
Krull dimension of R equals that of Rq. Hence the Krull dimension of G^ 
equals that of Gq. Since both G^ and Gq are connected and G^ C Gq, by 
the proposition in [T71 p. 25] we have G*^ = Gq. □ 

From Lemma [2.281 Gal(i2/A;) can be viewed as a subgroup of Gal(i?o/^'o)- In 
the following, we prove that Gal{R/k) is a normal subgroup of Gal(i?o/A;o)- 
Let J-'q and be the total ring of fractions of Rq and R respectively. Note 
that G8l{R/k) = Ga\{J^/k) and Gal(i?o/^"o) = Gal(J^o/fco)- Corollary [222] 
allows us to assume that Tq C J^. 

Lemma 2.30 Let u £ k be of degree m over /cq and let u = ui,U2, . . . , Um 
he its conjugates. Then there exist M E GLm(/co) and N G Qlm{kQ) such that 



Z 



(I 1 ... 1 \ 

Ui U2 ... Um 

\ ur' . . . ^ / 



satisfies 



a{Z) = MZ and 5{Z) = NZ. 



Proof. We claim that any field automorphism r of A; over kQ is a a5- 
field automorphism. Indeed, since k is an algebraic extension of kQ so any 
automorphism r is automatically a (5-field automorphism. One sees that r 
is a <T-field automorphism by noting that for any / G C(i)(x), r acts on 
the coefficients of powers of x while a acts only on x. For any g in the 
automorphism group of k over kQ, we have g{Z) = Z[g\ where [g] is a 
permutation matrix. Since g is also an automorphism of cr 5-fields, both M = 
a{Z)Z~^ and N = 5{Z)Z~^ are left invariant by g and therefore must have 
entries in fco. □ 

We now proceed to prove the main result of this section. 
Proposition 2.31 Gal{R/k) is a normal subgroup o/ Gal(i?o/A;o)- 
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Proof. Consider the following diagram 

T 

I \ 

k 

\ I 

J'onk 

ko 

First, we claim that the map Gal(jF/fc) Gal(.Fo/J^o H k) that sends h £ 
Gal(J^/k) to its restriction h\jr^ on J^q is an isomorphism. Any automor- 
phism of over k is determined by its action on a fundamental matrix 
of {o"(y) = AY,5{Y) = BY} and its restriction on J^q is determined in 
the same way. This implies that the restricted map is injective. To see 
that it is surjective, note that its image is closed and has .T-q H A; as a fixed 
field. Therefore, the Galois theory implies that the restricted map must 
be Gal(J^o/-^onA:). 

We now show that ghg~^ is in Gal{J^o/J^o n k) for any h G Gal{J^o/J^o n k) 
and g G Gal(.7^o/^o)- It suffices to show that g leaves J-'q Ci k invariant, 
which will imply that ghg^^{u) = u for any u £ J-Qr\k and so ghg^^ £ 
Gal(^o/-^o n k). Now let u £ TqCi k be of degree m over /cq- From Lemma 
12.301 U = (1, n, n^, . . . , satisfies some difference-differential system 

over ko- Therefore the vector g{U) satisfies the same system and so must be 
a C-linear combination of the columns of Z. In particular, we have g{u) E k. 
Therefore g leaves J-oDk invariant. This completes the proof. □ 

Theorem 2.32 If {a{Y) = AY,6{Y) = BY} is irreducible over ko, then it 
is equivalent over ko := FoDk to the system 

a{Y) = dmg{Ai,A2,--- ,Ad)Y, d(y) = diag(Si, ^2, • • • ,Bd)Y 

where Ai E GL£(A:o),-Bi E Qle{ko) and i = ^. Moreover, the system {o-{Y) = 
AiY, d(Y) = BiY} is irreducible over k for i = 1, . . . ,d, and there exists gi E 
Gal(i?o/^o) such that gi{Ai) = Ai and gi{Bi) = Bi. 

Proof. By Proposition 12.311 Gal{R/k) is isomorphic to Gal(jFo/A;o) and 
then Gal(.Fo/A;o) is a normal subgroup of Gal(i2o/^o)- Let V be the solution 
space in i?^ of{a{Y) = AY,6{Y) = BY}. Then Cliff'ord's Theorem ^ p.25. 
Theorem 2.2] tells us that V can be decomposed into F = Fi © V2 © • • • ®Vd 
where the Vi are minimal Gal(.7ij/A:o)-invariant subspaces of V and, for 
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each i, there exists gi G Gal(i?o/^o) such that Vi = gi{Vi). Further- 
more, g G Gal(i?o/A;o) permutes the Vi. Let Zi he an n x £ matrix whose 
columns are the solutions in Vi of the original system. Then Zi has the full 
rank. Then for each i, the columns of Zi = gi{Zi) are the solutions in Vi of 
the original system and Z^ has the full rank too. Let Z = (Zi,--- ,Z(i). 
Then Z is a fundamental matrix of the original system. By Lemma 1 
in [11], there exists i x n matrix Pi of the rank i with entries from ko 
such that PiVi = for i = 2, • • • ,d. Since g G Gal(i?o/^o) permutes the Vi, 

{9i{V2),--- ,gi{Vd)} = {Vi,--- ,V.i,V+i,--- ,Vd}. 

Let Pi = gi{Pi). Then PiVj = for j / i. Therefore, P := (Pf , • • • , Pj)^ G 
Qln{ko) satisfies that 

PZ = diag{Ui,U2,--- ,Ud) (2) 

with Ui G qIi{J-'q). Moreover we have Ui = gi{Ui) for each i. We now prove 
that det(P) 7^ 0. Assume the contrary that det(P) = 0. Then w^P = 
for some nonzero w ^ k^. Therefore there exists Wi G k^ for 1 < i < d 
such that wJ^Pi + • • • + wJPd = 0. Since the Pi have full rank, there exists 
at least one i such that wfPi 7^ 0. Without loss of generality, assume 
that wJPi / 0. From wJPi = - {wJP2 + ■■■ + w^Pa) and PiZi = for i = 
2, • • • , d, we have wJPiZ = 0. Since det(Z) ^ 0, wJPi = 0, a contradiction. 
Therefore det(P) / 0. Let i = ^. From Q, {a{Y) = AY,6{Y) = BY} is 
equivalent over kQ to 

a{Y) = diag(yli ,Ad)Y, 6{Y) = diag(i?i ,Bd)Y 

where Ai G GL^(A;o) and Bi G Qhiko)- Furthermore, Ui is a fundamental 
matrix of the system = AiZ,6{Z) = BiZ}. Since Ui = gi{Ui), we 

have that Ai = gi{Ai) and Bi = gi{Bi) for each i. From the minimality 
of Vi, the system {a{Z) = AiZ,6{Z) = BiZ} is irreducible over k. □ 

Corollary 2.33 Let Ai and Bi be as in Theorem [2JE for i = 1,... ,d. Then 
for each i, the Galois group of {a{Z) = AiZ,6{Z) = BiZ} over k is solvable 
by finite if and only if the Galois group of {a{Y) = AY, 5{Y) = BY} over ko 
is solvable by finite. 

3 Systems 

Throughout this section, let ko be the field C(t, x) with an automorphism a : 
X 1-^ X + 1 and a derivation ^ = and let k the extension field C(t)(x). In 
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this section, we will first prove that if a system {cr{Y) = AY,d{Y) = BY} 
with A G GL„(/i;o) and B G 0^n(^o) is irreducible over ko and its Galois 
group over ko is solvable by finite then there exists £ G Z>o with i\n such 
that the solution space of {a^{Y) = A£Y,6{Y) = BY} has a basis each of 
whose members is the interlacing of hypergeometric solutions over k. We 
will then refine this result to show that {a^{Y) = AiY,6{Y) = BY} is 
equivalent over k to a special form. Based on this special form, we will 
describe a decision procedure to find its solutions. 



3.1 Systems with Liouvillian Sequences as Solutions 



By Theorem ESU if a system {a{Y) = AY, 6{Y) = BY} with A G GL„(A;o) 
and B G Qlniko) is irreducible over /cq, then it can be decomposed into 
factors that are irreducible over k and if the Galois group of the system 
over ko is solvable by finite then the Galois groups of these factors over k 
are also solvable by finite. Hence it is enough to consider factors of the 
original system over k. 

Proposition 3.1 Suppose that {(t{Y) = AY,d{Y) = BY} with A G GL^(fc) 
and B G Qli{k) is an irreducible system over k and that its Galois group 
over k is solvable by finite. Then the system is equivalent over k to 



a(Y) = AY, 5{Y) = BY 



where B G Qli{k) and 



A 



/o 1 ••• o\ 

1 ••• 

••• 1 

\a ••• 0/ 



G GLeik) 



with a £ k. 



Proof. The proof is similar to those of Lemma 4.1 and Theorem 5.1 in \X3\- □ 



Remark 3.2 From the proof of Lemma 4-1 in fl3^, we know that £ di- 
vides |Gal(7^/fc)/Gal(7^/A;)°| because {a{Y) = AY,d{Y) = BY} is irre- 
ducible over k. From the proof of Theorem 5.1 in llSij . Gal(7?./fc)'^ is di- 
agonalizable. 



As a consequence of Proposition 13. H we have the following 
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Corollary 3.3 // a system {a{Y) = AY,6{Y) = BY} with A G GLeik) 
and B G Qle{k) is irreducible over k and its Galois group over k is solvable 
by finite, then {a^{Y) = AeY, 6{Y) = BY} is equivalent over k to 

a\Y) = VY, 6{Y) = BY 

where B G gle{k) and V = diag(a, (T(a), • • • , (T^"^(a)) with a as indicated in 
Proposition \3.1\ . 

Next, we shall prove further that {a^iY) = A(Y,5{Y) = BY} is equivalent 
over A; to a system of diagonal form. Note that equivalent systems have the 
same Picard-Vessiot extension. 

Proposition 3.4 // a system {a{Y) = AY,5{Y) = BY} with A G Ghi{k) 
and B G Qh{k) is irreducible over k and its Galois group over k is solvable 
by finite, then {a^{Y) = A£Y,6{Y) = BY} is equivalent over k to 

a\Y) = VY, 6{Y) = diag(fei , • • • , bi)Y 

with T> as in Corollary \3.3\ and bi £ k for i = 

Proof. By Corollary [331 W{Y) = AeY,5(Y) = BY} is equivalent over k 
to {a^{Y) = VY,d{Y) = BY} where B G gkik) and 

P = diag(a, cr(a), • • • ,a^~^{a)) 

with a as in Proposition l3.ll Let TZ = eo7?.©ei7^©- • -©e^T^ be the decompo- 
sition of 7^. Then Gal{TZ/k)^ = Ga\{eoTZ/k) by Lemma[22]and Gal(eo7^//c) 
is diagonalizable by Remark 13.21 From Lemma 12.31 it follows that cqTZ is 
a a'"S-P\ extension of k for the system 

a''{Y) = a'^-\V)---a{V)VY, 5{Y) = BY. (3) 

Let V = aJ-^{V)---a(V)V = diag((Ji,--- ,d(>). By Lemma [221 eo7^ is a 
domain. As in the differential case, we can show that ([3|) is equivalent 
over k to the system 

a''{Y)= diag(ai , • • • , a,)y, 6{Y) = diag(6i , • • • , 6^ )y (4) 

where aj,6j G k. Then there exists G = {gij) G GL^(/c) such that 

a"(G)diag(oi,--- ,ae)=VG, 
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which imphes (T^{gij)aj = gijdj. Since det(G) / 0, there is a permuta- 
tion ii, • • • , i£ of {1, 2, • • • ,i} such that gii^g2i2 • • • gue / 0. Hence 

(T^igjij) _ 
dj = Oi . for j = 1, • • • , 

Let P be a multiplication of some permutation matrices such that 

P^Miag(ai, • • • ,a^)P = diag(aii, • • • ^^ij 

and let T = Pdiag(l/5ij^ , • • • , l/gui)- Under the transformation Y TY, 
the system (j4]), and therefore ([3]), is equivalent over /c to 

a"(y)=W, (5(y) =diag(6i,--- 

where bi £ k for ? = 1,...,^. Proposition 12.261 implies that the sys- 
tem {a^(Y) = VY, d{Y) = BY} is equivalent over k to 

a\Y) = VY, 6{Y) = diag(6i, • • • , be)Y. 

This concludes the proposition. □ 

Theorem 12.321 together with Proposition 13.41 leads to the following 

Proposition 3.5 If the system {cr{Y) = AY, 6{Y) = BY} with A £ GL„(A;o) 
and B G 0/n(^o) is irreducible over ko and its Galois group over ko is solv- 
able by finite, then there exists i G Z>o with i\n such that the solution space 
of {a^(Y) = AiY,S(Y) = BY} has a basis consisting of the interlacing of 
hypergeometric solutions over k. 

Proof. By Theorem 12.321 and Proposition 13.41 there is ^ G Z>o with i\n 
such that {cr^(Y) = AiY,6{Y) = BY} is equivalent over A; to a system of 
diagonal form. Since the solution space of the latter system has a basis 
consisting of the interlacing of hypergeometric solutions over k, so does the 
solution space of {a^Y) = A^Y, 5{Y) = BY}. □ 

Corollary 3.6 Let C be the ring of liouvillian sequences over k. Assume 
that {a{Y) = AY,5{Y) = BY} with A G GL„(A;o) and B £ 0/„(A;o) is 
irreducible over ko and its Galois group over ko is solvable by finite. Then 
the solution space of the system has a basis with entries in C. 

Proof. Proposition 13.51 implies that there is ^ G Z>o such that the solution 
space of {a^{Y) = A£Y,S{Y) = BY} has a basis with entries in C. The 
corollary then follows from Proposition 12.241 □ 
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Let us turn to a general case where a difference-differential system may be 
reducible over the base field. If the Galois group over the base field of the 
given system is solvable by finite, then the Galois group over the base field 
of each factor is of the same type. The method in [13] together with the 
results in [7] implies the following 

Proposition 3.7 If the Galois group for {a {Y) = AY,5{Y) = BY} with A £ 
GL„(fco) and B G 0/n(^o) solvable by finite, then the solution space of the 
system has a basis with entries in C. 

Proof. By induction, we only need to prove the proposition for the case 
where the given system has two irreducible factors over ko. In this case, the 
given system is equivalent over /cq to 



where the systems {(7{Y) = AiY,6(Y) = BiY} for i = 1,2 are both irre- 
ducible over kf). Let di be the order of Ai for i = 1,2. By Corollarv l3.6^ each 
system {cr{Y) = AiY,6{Y) = BiY} has a fundamental matrix Ui £ Ghd^iC). 
From the proof of Theorem 3 in , Proposition 12.111 and Remark 12.121 it 
follows that the original system has a uJ-PV extension TZ of ko which con- 
tains entries of the C/j's and, moreover, has a fundamental matrix over 5^- 
of the form 



So a{V) = AiUi + A2V. Let V = U2W . Then a{W) = W + a{U2)-^AiUi. 
Since Ui G GLd.(-C), the entries of W are in C and so are the entries of V . □ 

3.2 Normal Forms 

Assume that a system {a{Y) = AY, 5{Y) = BY} where A G GL„(A:o) 
and B G 0^n(^o) is irreducible over k^ and that its Galois group over k^ 
is solvable by finite. Theorem 12.321 and Proposition 13.41 imply that there 
exists I G Z>o with i\n such that {a^Y) = A(Y,5{Y) = BY} is equivalent 
over A; to a system of diagonal form. In this section, we will show further 
that the above system is equivalent over /c to a more special form when the 
order n of the original system is prime. 
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3.2.1 Normal Forms for General Systems 

Let us first review some notions and properties concerning rational solutions 
of difference equations. 

Definition 3.8 (cf. JJM Definition 6.1]) Let f = ^ with P,Q e C(t)[x] 
and gcd(P, Q) = 1. 

(1) The dispersion of Q, denoted by disp{Q) is 



max{j G Z>o|(5(a) = Q{a + j) = for some a G C(t)}. 

(2) The polar dispersion of f is the dispersion of Q and denoted pdisp(/). 

(3) / is said to be standard with respect to a"^ , with m G Z>o, i/disp(P • 
Q) < m. 

As in [12], we have the following 



Lemma 3.9 Assume that f £ k \ {0}, a G C{t) \ {0} and m G Z>o. 

(1) There exist f,g £ k \ {0} such that f = ^—^f where f is standard 
with respect to a"^. 

(2) /// has a pole, then pdisp((T™(/) — af) > m. 
Proof. The proof is similar to that of Lemma 6.2 in |12j . □ 



Proposition 3.10 Let ^ a,b £ k satisfy a"^{b) — b= where m G Z>o. 
Then 

/ / ml3{t) 



where f G k,c,(3{t) G C{t), and a{x) G C(x) is standard with respect to a"^ . 

Proof. Let a = ^—j^a with a standard with respect to o"'" and b = 6— ^y^. 
Then a'^ib) — b = View a and b as rational functions in x. Then 

pdisp(^fc^) < m. If ^ C(t), then has a pole and so does b. 
By Lemma 13. 9| pdisp(cj'"(6) — b) > m, a contradiction. Hence = 
uj{t) G C(t), which means that a = a(x)e-^"'^*^'^*. Since a G A:, a is of 
the form a{x)P{t) where a{x) G C(x) and (3{t) G C{t). Then b G C(t)[x]. 
Suppose that b = Cnx'^ + Cn-ix^~^ + • • • + cq where Cj G C(t) and / 0. 
Then 



So n = 1 and 6 = ^^^^x + cq. □ 
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Theorem 3.11 If {a{Y) = AY,6{Y) = BY} with A G GL„(A;o) and B G 
0^n(^o) is irreducible over and its Galois group over Uq is solvable by 
finite, then there exists I £ Z>o with £\n such that the system 

a'{Y) = A,Y, 5{Y) = BY 

is equivalent over k to 

' cT^(y) = diag(A(x)/3i(t),A(x)/52(t),-- - , A(x)/3^(t))y, 



m = diag mmxi, +c,,.--, + c^)) y 



(5) 

where A(x) = diag(a(x), • • • , a{x + £ — 1)), Ci = diag(ci, • • • , q) and mi=n. 
Moreover, a{x) G C{x) is standard with respect to a^, l3i{t),Ci G C(t), 
and there exists gi in the Galois group of the original system over ko such 
that I3i{t) = gi{l3i{t)) and d = gi{Ci). 

Proof. By Theorem 12. 32^ it suffices to prove the theorem for a factor of 
the given system. Let {criY) = AY, ^{Y) = BY} be such a factor with A G 
GU{k) and B G Qk{k). By Proposition [331 W^iY) = AiY,5{Y) = BY] is 
equivalent over k to 

cj\Y)=VY, <5(y) =diag(5i,--- 

where P is as in Corollary 13.31 and bi £ k for i = Since and 6 

commute, we have (7^{bi) — bi = By Proposition 13.101 we have 

a = — —a{x)Pi{t) and bi = —— + x + ci 

f J >^Pi\t} 



where a{x) G C{x) is standard with respect to a , ci,(3i{t) G C{t) and f £ k. 
Then for i = 1, • • • ,i, 

^ ^ = — i~i(f\ a{x+i-l)[3i{t) and bi = — + x+q. 

Let F = diag(/, cj(/), • • • ,a^~^{f)). Then the system 

a\Y)=VY, 5{Y) = dmg{bi,--- ,b,)Y 
is equivalent over k to 

a'{Y) = A{x)Pi{t)Y, 6{Y) = (^J§^xl, + Ci] Y 



under the transformation Y — > FY where A(x) = diag(a(x), • • • , a{x+i—l)) 
and Ci = diag(ci, • • • , q). □ 
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3.2.2 Normal Forms for Systems of Prime Order 



If a difference-differential system is of prime order n, then the integer £ in 
Theorem 13 ■ 11 1 eg uals either 1 or n. For the case where the system is reducible 
over k, we can refine Theorem 13.111 further in the following 

Proposition 3.12 Assume that n is a prime number. Suppose that the 
system {a(Y) = AY,6{Y) = BY} with A G GL„(/co) and B G 0/n(A:o) 
is irreducible over and reducible over k and that its Galois group of is 
solvable by finite. Then the system is equivalent over k to 



'a(y) = a(x)diag(/3i(t),/J2(t),--- , Pn{t))Y, 



5(y) = diag(^fflx + ci, 



iiS^x + Cn))Y 



(6) 



where a{x) G C{x) is standard with respect to cr, f3i{t) = gi{Pi{t)) G C{t) 
and Ci = gi{ci) G C(t) for some gi in the Galois group of the original system 
over kQ. 

Before discussing the other case where a difference-differential system is 
irreducible over k, let us look at the following 

Lemma 3.13 Assume that (t{Y) = AY with A G GL„(A:o) is equivalent 
over k to (j{Y) = AY where 



( 



A 










V/3(t)a(x) 



1 

1 







0\ 


1 

0/ 



with a{x) G C{x) and (3{t) G C{t). Then f3{t) G C(t). 
Proof. There exists G G GLn{k) such that a{G)AG-'^ = A. Then 
det(a(G)) det(l) det(G-i) = det(^). 



Since det(cj(G)) = a(det(G)) and det(G-i) 



1 



dct(G) 



we have 



(-ir-i/3(tHx)^^^^ = det(A). 



det(G) 

Expand the rational functions in x in the above equation as series at x = oo. 

(det(G) ' 
det{G) 



Since = ^ + iQ where Q G C(t)[[i]], one sees that (5{t) G C(t). □ 
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Proposition 3.14 Let A,A e GL„(A;o). If a{Y) = AY and a{Y) 
are equivalent over k then they are equivalent over /cq- 



AY 



AG. 



Proof. Suppose that there exists G G GL„(A;) such that a{G)A 
Then there exists 7(t) E C{t) such that G G GL„(/co(7(i)))- Let m = 
[ko{'y{t)) : ko]. Since l,^(t),--- ,7(t)"*""'^ is a basis of ko{'y{t)) over ko, we 
can write 

G = Go + Gi7(t) + • • • + G^_i7(t)™-i 

where Gj G Qlniko)- From cr(G)A = AG, it follows that cr{Gi)A = AGi 
for i = 0, ••• ,m — 1. Let A be a parameter satisfying o"(A) = A and 
let H{X) = T.T=o^X'Gi. Therefore, a{H{X))A = AH{\). Since det(G) = 
det(i?(7(i))) 7^ 0, det(-ff(A)) is a nonzero polynomial with coefficients in k^. 
Hence there exists c G C(t) such that det(i?(c)) / 0. So c7(i/(c))A = AH{c) 
and i?(c) G GL„(/co). □ 

We now turn to the case where a difference-differential system over k^ is 
irreducible over k. 

Proposition 3.15 Suppose that{a{Y) = AY,5{Y) = BY] with A G GL„(A;o) 
and B G 0/n(^o) irreducible over k and that its Galois group over ko is 
solvable by finite. Then the system is equivalent over ko to 



a{Y) = AY, 5{Y) = BY 



where B G Qln{ko) and 
( 

A = 







0\ 


1 

0/ 



G GL„(A;o) 



\(3{t)a{x) • 
with a{x) G C(x) standard with respect to a" and (3{t) G C(t). More- 

Proof. By Proposition 13. H the given system is equivalent over k to the 
system {a{Y) = AY, 6{Y) = BY} where B G and 



/O 1 
1 



A 



0\ 





\a 



1 

0/ 



G GLn{k) 
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for some a £ k. Since a and 6 commute, we have a{B)A 
Let 8 = (bij) 

nxn where bij G k. Then 
6{a) - 



6{A) + AB. 



a{b„ 



Hence cr"'(6ii^ 



a 

5{a) 



•J22 



By Proposition 13. lOl we have 



/ 



a{x)(5{t) 



with f £ k, a{x) G C{x) standard with respect to a" and (3{t) G C(t). 
Then {a{Y) = AY, 6{Y) = BY} is equivalent over k to {a{Y) = AY, 5{Y)_ = 
BY} under the transformation Y - 
Qln{ko) and 

/ 





diag(/, a{f), ■■■ , where B G 



A 







1 

1 









1 

0/ 



with a{x) G C{x) and (3{t) G C(i). By Lemma 13.131 and Proposition 13.141 
the original system is equivalent over k^ to {(j{Y) = AY,5(Y) = BY} 

for some b G C{x) and 



with f3{t) G C{t). Assume that 



let n = • • • a{b)b. We have 



a{x) 



u(x) 



thus a{x) 



cu{x) for 



some constant c with respect to a. Therefore c G C since a{x) and u{x) are 
both in C(x). Let P G GL„(C(t)) be such that 

/ 1 ••• 0\ 
1 ••• 



P 



-1 





\cl3{t) 



1 

0/ 



P = diag(/3i(t),/32 (*),••• ,/3n(t)) 



where the Pi{tys are the roots of Y"- — c(3{t). Let 

F = diag(l, b, a{b)b, ■■■ , a^'-^ib) ■ ■ ■ a{b)b)P 

and B = F^^BF - F^M(F). Then {a{Y) = AY, biJ) = BY} is equivalent 
over A; to 

cT(y) = 6- diag(/?i(t),/?2 (*),••• Jn{t))Y, 5{Y) = BY 

under the transformation Y — > FY. Assume B = {bij)nxn- Since a and 6 
commute, 

0, 
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for all i,j with 1 < i ^ j < n. Hence bij = if i 7^ j. In other word, B is of 
diagonal form. This contradicts to the irreducibility over k of the original 
system. □ 

Lemma 3.16 Let a G A:o \ {0}, n be a positive integer and m > be the 

least integer such that = for some b ^k^. Then m\n. 

Proof. Suppose that = with 6 G /cq. Then for each i > 0, 

— = — ^— with a £ ko. 

a Q 

Let n = iim + £2 where < £2 < m — 1. Then 

(7^2 (a) a^i™+^2(a) a^^{a) cj"(a) c 



a a cr^i'»+^2(a) a cr"(c) 

for some c G k^. Hence ^2 = and so m\n. □ 

Proposition 3.17 Assume that n is a prime number, the system 

a{Y) = AY, 6{Y) = BY 

with A G GL„(A:o) and B G 0/n(^o) is irreducible over k and its Galois group 
over /cq is solvable by finite. Then {a^ {Y) = AnY,5(Y) = BY} is equivalent 
over ko to 

'o-"(y) = ;g(t)diag(a(a;),--- ,a{x + n-l))Y 
S{Y) = (^^xln + dmg{bi,--- ,bn))Y 

where a{x) and P{t) are as in Proposition \ 3. 1 5\ and hi G C(i) for i = 1, . . . ,n. 

Proof. By Proposition Eig {cr"(y) = AnY,6{Y) = BY} is equivalent 
over ko to the system 

(T"(y) = P{t) • diag(a(x), • • • ,a{x + n-l))Y, 5{Y) = BY 

with a{x) and (3{t) as in Proposition 13.151 and B G Qln{ko). Let B={bij)nxn- 
From cr"5 = , we have 



m 

Hence ba = with hi G C(i). Note that n is prime and ~^^^y^7^~~^ 

for any b G C(x). Then by Lemma [3. 161 "^^^^'^ 7^ ^ < ^ < — 1 

and b G C(x). Hence bij = for i 7^ j. This concludes the proposition. □ 
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3.3 A Decision Procedure for Systems of Prime Order 



Consider a system {cr{Y) = AY,6{Y) = BY} over ko. Assume that the 
order n is prime, the system is irreducible over fco and its Galois group 
is solvable by finite (or, equivalently, the system has liouvillian solutions). 
By Proposition 13.121 and Proposition 13.171 either the original system has 
hypergeometric solutions over k or the system {a^(Y) = AY,6{Y) = BY} 
has solutions which are the interlacing of hypergeometric solutions over /co- 
in this section, we will give a decision procedure to find solutions of systems 
of both forms when the order n is prime. Our procedure relies on the 
following three facts in the ordinary cases: 

(Ai) we can compute all rational solutions in k"" of an ordinary difference 



equation a(Y) = AY where A G GL„(A:); ( P El El [H] ) ; 

(^2) we can compute all hypergeometric solutions over C(x) of an ordinary 
difference equation a{Y) = AY where A e GL„(C(x)) ( [IS [2ll [22l [231 

mm): 

(^3) we can compute all hyperexponential solutions over C(t) of an ordinary 
differential equation b{Y) = BY where B G GL„(C(t)) ([UlEIlEa 



In the following subsections, we will reduce the problem of finding solutions 
of {cj(y) = AY, bij) = BY} or of {cj"(y) = AY, 5{Y) = BY} to that in the 
ordinary cases as indicated above. We have two case distinctions according 
to the reducibihty of {a{Y) = AY, 5{Y) = BY} over k. 

3.3.1 The Decision Procedure for the Reducible Case 

Assume that {criY) = AY, 5{Y) = BY} is reducible over k. Proposition [3]T2] 
implies that this system has hypergeometric solutions of the form Wihi 
for i = 1, . . . , n, where Wi S A;" and hi satisfies 



with a{x) G C(x) standard with respect to a, Pi{t) = gi{l3i{t)) G C{t) 
and Cj = gi{ci) G C(t) for some gi in the Galois group of the original system 
over fco. Substituting each Wih into the original system, we get 



[271 US]). 




T{Wi) 
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So, to compute hypergeometric solutions of {(t{Y) = AY, 6{Y) = BY} it 
suffices to find a{x), Pi{t), Ci and Wi satisfying d?]). 

Remark 3.18 The equalities ([7]) still hold when replacing a{x) by 
and Wi by ^ for g G C(x). So in the sequel, we will compute a suit- 
able instead ofa{x). 



Computing a{x): By Proposition 13.121 there exists G G GL„(A;) such that 
a(det(G)) 



det(G) 



-a{xTYlPi{t) = det{A). 



Without loss of generality, we assume that the numerator and denominator 
of a{x) are monic. Expanding the functions in the above equality as series 
at X = oo, one can compute HILi Pii^) from the series expansion of det(74) 
at X = oo. Let a = j^^^^Pit) • Rewrite d = where b,d £ ko and a is 

standard with respect to a. Then 



9 



a(x)" for some g € ko. (8) 



From Proposition 13.121 a{x) is standard with respect to cr and so is a(x)". 
Proposition 13.20] below shows that S C(x) and thus d G C(x). More- 
over, o has the form (-^a(x)) for some g £ C(x). To prove Proposi- 



g 

tion 13.201 let us introduce a notation used in [251 Section 2.1]. 



Definition 3.19 A divisor D on F^{C{t)) is defined to be a finite formal ex- 
pression Uplp] with p £ P^(C(t)) and rip £ Z. The support of a divisor D, 
denoted supp(L'), is the finite set of all p with Up / 0. Let p S supp(Z)). 
The Z-orbit E of p in supp(L') is defined to be 

E{p, supp(L>)) = {p + i\i £ Z and p + i £ supp(L')}. 

As usual, the divisor div(/) of a rational function f £ k \ {0} is given 
by div(/) = ordp(/)[p], where ordp(/) denotes the order of / at the 
point p. It is clear that dw{fg) = div(/) + div(g'). Moreover, if p is 
in supp(div(/)) but not in supp(div(/g)), thenp £ supp(div(/))nsupp(div(5')). 
By Definition 13.191 if / £ k \ {0} is standard with respect to a, then 
supp(div(/))) = {p} for each p £ supp(div(/)). 
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Proposition 3.20 Assume that f,gG ^\{0} and f is standard with respect 
to a. If o'{g)g~^ f is standard with respect to a, then 

o-(g) _ T-r {x + kj- Cj)^^ 
g ~ -l-i (x - Ci)-^^ 

with ki £ Z, nii £ Z>o, Q G C(t) and disp{Ylj^{x — Cj)) = 0. Moreover, for 
each i, either ordaif) = rui or orda-kiif) = -m- 

Proof. Let H = a{g)g~~^f, Si = supp(div(/)), S2=supp{dW{a{g)g-^)) 
and S3 = supp(div(i?)). By Lemma 2.1 in [25], 

ordg ( = for each p € 
qeEipA) V f / 

Then \E{p, 5*2)! > 2 for each p £ supp(S'2). Since H and / are standard, 

\E{p, S2) n < 1 and \E{p, ^2) n Si| < 1 

thus \E{p, S2) n {Si U Ss)! < 2. From ^2 ^ U S3, we have \E{p, 82)] < 2. 
Hence for each p £ S2, 

\E{p, 52)1 = 2, \E{p, S2) n 5i| = 1 and \E{p, S2) n ^3! = 1. 

From \E{p, S2)\ = 2 and \E{p, 52)n53| = 1, either ordp{a{g)g^^) = — ordp(/) 
or ordp+jo ((7(5)5'^^) = -ordp+jo(/) with p + jo ^ E{p,S2). The proposition 
holds. □ 

Let g be as in ([8]). Since a{x) £ ^{x), g can be chosen in C(x) according to 
Proposition 13.201 Then a £ C(x). Moreover, 

(y{g) _jT {x + ki- q)""' 

g (x — Cj)'"' 

where mj has the form ifiin for some mj £ Z>o since mj is either ordc^ (a(x)") 
or -orde,-fc, {a{xY). Let 5 = Jli EliLV + J " CiT' ■ Then 

Note that the numerator and the denominator of a{x) are monic, so we can 
compute ^^a{x) from a. 
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Example 3.21 Consider the integrable system 

a{Y) = AY, 6{Y) = BY 

where 

{x'^+\){f^t^-x'^-x) t{x'^^\)(f~x) 



B 



We have 



2xt^-t'^-t'^+t^+3t^+2t+x^t+t'^x+x 

-f^x^+t'^+2t'*+t^-x^+2t^x+a: -x^t-t^x-x+t^+2t'^-xt+t'^+t^ 
(f^-x){f^+l) {f^-x){f^+l) 



— j; — 1 

T/iiis /3i(t)/32(i) = -(t^ + 1). Letd = and wriie 

x2 + 1)^ 



^ ' 2 I i\2 



t2 - (x + 1) 

T/ien a(x) = + 1. 



Computing We first prove the following 

Lemma 3.22 Either f3i{t) = ■■■ = G C(t) or Pi{t),--- , Pnit) are 

the conjugate roots of an irreducible polynomial of degree n with coefficients 
in C{t). 

Proof. Let Ro be a a5-P\ extension of ko for {a{Y) = AY,5{Y)=BY} 
and P = nr=i("'^ ~ Piii))- From the proof of Theorem 12.321 one sees 
that Gal{Ro/ko) permutes the f3i{t). Furthermore, the orbits of the Pi{t) 
under this group action all have the same size. Therefore, P is a polynomial 
with coefficients in C(t). Since n is prime, either P is irreducible or all the 
factors of P in C(t)[X] are of degree one. This concludes the lemma. □ 

The following two notions can be found in [U [5] . 

Definition 3.23 Let H = {hij)nxn £ GL„(A;o)- The order of H at oo is 
defined as 

ordool-f^) = min{ordoo(/ijj)} 
where ordoo(^ij) is the order of hij at oo. 
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We rewrite H into the form 



H={-] IHo + Hi- + 



X 



where Hi G g/„(C(t)) and Hq / 0. 
Definition 3.24 The rational number 

m{H) = -oid^{H) + 



rank(i7( 



n 



is called the first Moser order of H . And 

fi{H) = mhi{m{a{G)HG~^)\G G GL„(A;)} 

is called the Moser invariant of H . A matrix H is said to he irreducible 
if m{H) = fJ-{H), otherwise it is called reducible. 

Given H £ GL„(A;o)i one can use the algorithm in [H [5] to compute G G 
GL„(fco) such that H := a{G)HG~^ is irreducible. So we can assume 
that is irreducible where A and a{x) are as in ([7]). Under this as- 
sumption, we will show that 

^0+^1- + 



a[x) X 

with Ai G gZ„(C(t)) for each i and that ah the Pi{tys are eigenvalues of Aq. 
The following lemma can be deduced from the results of Barkatou in [5]. 
We will present a self contained proof due to Reinhart Shaefke. 

Lemma 3.25 Let G G GL„(A;) and assume thatoidooi(7iG^^)G) =0. Then 
all the eigenvalues of a{G~^)G\x=oa are 1. 



Proof, hei H = a{G"^)G. Then H = Ho + Hi^ + ■ ■ ■ with Hi £ gln{C{t)) 
and Ho ^ 0. We now show that Hq — In is nilpotent. For a positive inte- 
ger m, consider a map Lm ■ Qln{k) — > Qln{k) given hy U ^ a{U)—a"^{H)U 
for any U G Qln{k). Set Pm = Lm o L^-i o • • • o Lo(/n) where o denotes 
the composition of maps. Then PmU=oo = {Ln — Hq)^~^^. On the other 
hand, Lm(c7™(G-^)y) = o-™+i(G-^)A(y) where A = cr - 1 is a difference 
operator and V G 0ln{k). Hence P„, = a'^+^{G-^)A"'+^{G). Note that 
when m increases, ordoo(A™+^(G)) increases but ordoo(c™'''^(G'~"^)) is in- 
variant. Then for a sufficiently large m, Pm\x=oo = 0. This concludes the 
lemma. □ 
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Now we can prove the following 

A 
a{x) 



Proposition 3.26 ordoo ( -r-r ) = and (3i{t) , . . . , (3n{t) are eigenvalues 



A I 

a{x) \x=oo' 



of 

Proof. By Proposition 13.121 there exists G G GL„(A;) such that 

a(G)^G-i = diag(/3i(t),--- 

a[x) 

This implies that ordoo (det (^)) = and m = < 1. By 

the property of orders, 

Since m (^-^^^ ^ 1) ordoo = by the definition of the first Moser 

orders. Therefore, 

^0+^1- + ••• 



a{x) X 

where Ai £ 0/„(C(t)) and Aq / 0. From i^, a{Y) = ^(^^y^^^t^ Y has a 
rational solution Wi in /c". Suppose that 

/-.N ordoo (W'i) / I 



where Wij G C(t) and Wio + 0. Then Wiq = -^W-ia. Since Wio / 0, 

we have det i^In — = 0- Hence all the Pi{t) are the eigenvalues of ^o- 

If the (3i{t) are the conjugate roots of some irreducible polynomial with de- 
gree n, then they are clearly eigenvalues of Aq. Thus by Lemma [3.22l we only 
need to consider the case f3i{t) = ■ ■ ■ = (3n{t) G C(x). In this case, = 

(3i{t)a{G~^)G. Since ordoo (^) = 0, we have ordoo (ct(G^I)G) = 0. By 

Lemma [3.251 all the eigenvalues of (j{G~^)G\x=oa equal 1. Hence all the 
eigenvalues of equal Pi{t). □ 



Example 3.27 (Continued) Let A = ^^^^ . From the process in we 
can find an irreducible matrix A which is equivalent to A where 



A 



{x+l)t(t^+l-x) x+1 
{t'-^-x-l)x f-'-x-l 

{t'-^-x-l)x t'^-x-1 
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Write A = Aq + Ai^ -\ where 

io=("/ J) and A^=(^, 
The eigenvalues of Aq are 

Computing Cj and Wj: Let A(t) = diag(/3i(f), . . . ,l3n{t)). From (^i) we 
can find a matrix G £ GL„(A:) such that a{G)a{x)A{t) = AG. Let B = 
G-^BG-G-^6{G). Then Be Qln{k) and the system {a{Y)= AY, 5{Y)=BY} 
is equivalent over k to 

a{Y) = a{x)A{t)Y, 6{Y) = BY. (9) 

Note that G may not be the required transformation matrix in Proposi- 
tion 13.121 so B may not be of diagonal form. Since a5 = 6a, the same 
argument as in the proof of Proposition 13.151 implies the following conclu- 
sions: 

(i) If (3i{t) / (3j{t) for all i,j with I < i ^ j < n, then 



with Cj S C{t); 

[ii) If I3i{t) = ■■■ = (3n{t) G C(t), then by Proposition EUl G can be 
chosen in GL„(fco). Thus 

with B e QlniC{t)). 

In the case (i), we obtain the Cj's, and the Wi's are just the columns of G. 
For the case (ii), since ([9]) is equivalent over k to there exists G G GL„(A;) 
such that 



(TiG) = G, 

6{G)+G('-f{fxIn + diag{c,,--- ,cn)) = [b + '-f^xl^r) G. 



Hence G £ GL„(C(t)) and 5{Y) = BY is equivalent over C{t) to 

5{Y) = diag(ci,--- ,c„)y. 

Solving the system 6{Y) = BY by (^43), we get the Cj's and the Wi are just 
the columns of GG. 
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Example 3.28 (Continued) Let K{t) = dmg{\/t'^ + 1, -VWTl). From (Ai), 
we can obtain G G GL2(A;) such that a{G){x'^ + l)A(t) = AG where 



G 



-x+Wf^+l _x±tVE±l I 
2{f^~x) 2{t'^-x) / 



Then 



B = G-^BG-G-^5{G) 



xt 



^-V^^+ij- 

Hence ci = + 1 + 1; C2 = + 1 and Wi is the i-th column of G 

for i = 1,2. Furthermore, a basis of the solution space is 



I tV¥+i~x I \ x+tyfWyL 

\ 2{t'^-x) I \ 2{x-f^) 

where h satisfies that cr{h) = {x^ + \)h and 5{h) = 0. 

3.3.2 The Decision Procedure for the Irreducible Case 

Assume that {a{Y) = AY,5{Y) = BY} with A G GL„(A;o) and B £ 0/„(/eo) 
is an irreducible system over k and its Galois group over ko is solvable by 
finite. By Proposition EUl the system {o-"(y) = AnY,6{Y) = BY} has 
solutions of the form Wj/ij for i = 1, • • • ,n, where Wi G k^ and hi satisfies 

a^ih,) = aix + i- l)Pit)hi, 5{hi) = + 

with a{x),(3{t) and bi as in Proposition 13.171 Substituting Y = Wihi 
into {(J''{Y) = AnY, 6{Y) = BY}, we have 

(10) 

To compute Wihi, it suffices to compute a{x), l3{t),Wi and bi which sat- 
isfy (jlOp . Without loss of generality, we assume that the numerator and de- 
nominator of a{x) are monic. By Proposition 13.151 there exists G G GL„(/eo) 
such that 

^^^(-ir-ia(x)/3(t)=det(^). 
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Expanding det{A) as a series in we get that (— is the leading 
coefficient of the series. Hence we can obtain (3{t) from det(^). In this 
case, we can not find a{x) by the method used in Section 13.3.11 However 
we can reduce this problem to working with difference equations over C{x). 
By Proposition 13.171 there exists G G GL„(A:o) (the same as that in Propo- 
sition I3.15P such that 

• diag(a(x), . . . , a(x + n - 1)) = -^G. 

f3{t) 

Assume that t = p \s not a pole of the entries of and such that 
det (^\t=j) / 0. Let ^ = AQ + {t-p)Ai + - ■ ■ where Ai G Qln{C{x)). We 



will show that a{x) can be found by examining the hyper geometric solutions 
of 0""(y) = AqY . This will follow from the next proposition. 

Proposition 3.29 Some factor ofa^iY) = AqY is equivalent overC{x) to 
some factor of a"'{Y) = diag(a(2;), a{x + 1), • • • , a{x + n — 1))^. 

Proof. Let G be as above and let ^'(x) = diag(a(3;), • • • ,a{x + n—l)). We 
may multiply G by a power of t —p and assume that G = Go + {t—p)Gi + - ■ ■ 
where Go / and Gi G Qln{C{x)). Then 

ct"(Go + {t- p)Gi + ■■■ )^{x) = (Ao + . . . )(Go + (t - p)Gi + •••). 

Therefore <t"(G'o)^(x) = AqGq. Let r = rank(G'o). Then r > because Gq/O. 
There exist P G GL„(C(x)) and Q which is a product of some permutation 
matrices such that 



G = PGoQ 





G21 G22 



where G22 G GL^(C(x)). Then 

a"(G)diag(a(x + fci), • • • , a(x + K)) = a''{P)AQP-^G (11) 
where ki, - ■ ■ ,kn are a permutation of {0, 1, • • • , n — 1}. Now let 



A = a^{P)AoP-' = { ] where A22 G b/.(C(x)), 

\A2l A22/ 

and D2 = diag(Q(x+/c„_r+i)5 • • • , a(x+/c„)). From (fTTj) . we have ^i2G22=0 
and a"(G22)-D2 = ^22622- Since G22 G GL,.(C(x)), we have Au = 0. There- 
fore ct'^{Z) = A22Z is a factor oia^{Y) = AqY, which is equivalent over C(x) 
to o-"(Z) = D2Z. □ 
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Remark 3.30 For almost all of p £ C, a'^{Y) = -^\t=pY is equivalent 
over C(x) to 

cr"(y) = diag{a{x),a{x + !),••• ,a{x + n- 1))Y. 
since G\t=p is invertible. 

The same argument as in Remark 13.181 imphes that it is enough to com- 
pute ^9)ot{x) £qj, some suitable g G C(x) instead of a{x). We can use 

Proposition 13.291 to find °" (g)^(^+^) -v^rith k £ Z and g G C(x) as fohows. 

From Theorem 3 in [7], if (-^i,-"" i-^r)^ is a solution of cr"'(Z) = A22Z, 
then (0, • • • , 0, zi, • • • , ZrY^ is a solution of C7''(y)_ = AY. So <T"(y)_= AqY 
has at least r solutions Wihi, . . . , Wrhr, where Wi G C(x)"' and /ij satis- 
fies (T"(/ij) = a{x + kn-r+i)hi. By (^2)1 we can find all hyper geometric solu- 
tions oia{Z) = AQ{nx)Z where AQ{nx) means replacing x by nx in Aq. Then 
by interlacing, we can find all solutions of a^iY) = AqY of the form Wjhj 
where Wj G C(x)" and hj satisfies a^{hj) = djhj for some dj G C(x). Then 
there exists hj^ such that hj^ = ghi for some g G C(x) and 

0(2; + = : = -^a{x + kn-r+l). 

hj, 9 

After finding a{x + kn-r+i), we can compute a matrix G G GL,„(A:o) in a 
finite number of steps by (^1), such that 

a''{G-^)AnG = /3(t)diag(a(x), ■■■ ,a{x + n- 1)). 

Let B = G^^BG — G~^5{G). Then we get a new system 

fj'^(y) = /?(t)diag(a(x), ■■■ ,a{x + n- 1))Y, 6{Y) = BY 

which is equivalent to the original one under the transformation Y — > G~^Y. 
Since <t" and 6 commute and "^^^-^^^ 7^ ^ ^ C!(a;), the same 

argument as in the proof of Proposition 13.171 implies that B is of diagonal 
form, that is 

B = diag ^^TV^ + &!,•••, ^^TT^ + 
We then get the bi, and the Wi are just the i-th columns of G. 
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Example 3.31 Consider an integrahle system: 

a(Y) = AY, S{Y) = BY 

where 



/ xH'^+2xH^+xt^~x~l t^{tx^+2tx^+tx^+l) t{t-x~l) \ 



A 



\ 



W+x+l 

£x^+^x+x+l 
t{fi+x+l) 



B 



( t'''+t'^x+x^+t*x-t^ 
' t{fi+x) 

t'-'+x 



t^+x 



W+x+T 
T^Tx+i 

t'^+x+1 

xj-t'^+t^-l) 
fi+x 



2t'^x+x^+t^-t'^ 
t(W+x) 

x(-t^+f^-l) 
t(W+x) 



t'^+x+1 

t(l+t) 

W+x+1 

t-x-1 
W+x+l / 

xt^{-l+t) 

_t_H-i+ti 

t'^+x 



x^+xt^+t'^x+f-x-t" 



6 



i(W+xY 



We have 



+ X + 1) 



-By ^/le Galois group over ko of the given system is solvable by finite, 

then this system have no hypergeometric solutions over k. Therefore we 
consider the system 

a^{Y)=AsY, 6{Y) = BY 



where 



A. 



/ t^(t'^x^+t^x+21x+x^ +8x^+18) t'^{x+l){5x+6) 
I t'^+x+3 t'^+x+3 



V 



2t^2x+3) 
t'^+x+3 

2t*{2x+3) 
t'^+x+3 



{x+l)t^{x'^+t'^x+2t'^) 
t'^+x+3 

(x+l)t^{5x+6) 



2t4(x+2)(x+3) \ 
P+x+3 > 

2t^{x+2) 
t'^+x+3 

(x+2)(x+3)t3 (x+l+t^) 
fi+x+3 ' 



We can compute (3{t) = t'^ from det(^). Let A= Then 



t=o 



/{x + 2){x + 3) 

Q {X+1)X^ Q 

V ^^^4S^ ix + l){x + 2)) 
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By {A2), all hypergeometric solutions of a^{Y) = A\t=oY are 



^0^ 



9^r(^)r(^) , 9fr(^)r(^) , 9lr(|)r(^) 



.0, 



,1, 



where T{x) satisfies T{x + 1) = xT{x). By (Ai), we can compute a rational 
solution of a^(Y) = ^(^^"^i^i Y- Moreover, we can compute a matrix G G 
GL3(C(x,t)) such that 

o-^(G)diag(x(x + {x + l)(x + 2)t\ {x + 2){x + 3)t^) = A3G 

where 

G = 



— X 

t t 

X t 



'F+x t(F+x) F+x / 

Let B = G-^BG - G~^5{G). Then B = diag (f + f + + t^) ■ Hence 
a basis of solution space of {a^{Y) = A3Y,6{Y) = BY} is 



F+S 



Fi(x):=9tr(|)r(^)t^e^ ^ 



V2ix) := 9f r (2=±i) r (^) t^e- 



Viix) :=9tr(^)r(^)i-eT 



'F+x/ 



— —JL— 

F+x 

F+x 

X 

^t(F+x)^ 



F+x 
t 

F+x 



t^+x 



Clearly, Vi{l) 7^ /or i = 1,2,3, and j4(j) and -B(j) are well defined 
and det(A(j)) 7^ for j > 1. 5?/ the results in Section lK^ we get a basis of 
the solution space of the original system: 



/in t it^ _ 

/v^'F+T'F+2' F+3 



Wi = mTa)Ta)te^ 



1 



2f* 2t* 



'^^'F+T' F+2'F+3 



1 'It'^ 

F+T' F+2' F+3' 



•••A 
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and 



/ V'J) 75+2' FT3' 



W2 = 93r(f )r(i)te 



W3 = mT{l)T{^)te- 



+2' t2+3' 
V (0 1 ^ 6*4 , 



(0, 



+1' 7^+2'F+3' 



A'"oie i/iai aZZ ifte W,' are liouvillian. 



•) 

•v 

) \ 

) 

•••)/ 



3.4 Summary 

Consider two systems 



a{Y) = AY, 6iY) = BY 



and 



AnY, 



6(Y) = BY 



(12) 



(13) 



a-{Y) 

where A £ GL,j(A;o)) B £ Qlniko) and n is a prime number. Assume that (jl2p 
is irreducible over ko . From the results in Sections 13.3.11 and 13.3.21 if (fT2|) 
has a houvillian solution over k, then either the solution space of (|12p has 
a basis consisting of hypergeometric solutions over k or the solution space 
of (jlSp has a basis each of whose members is the interlacing of hypergeo- 
metric vectors over ko. Let us summarize the previous decision procedure 
as follows. 



Decision Procedure 1 Compute a fundamental matrix of (|12p whose en- 
tries are hypergeometric over k if it exists. 

(a) Write det(A) = where g,a £ ko and a is standard with respect 
to a. If a ^ a(xf"/3(t) for any a{x) G C(x) and P{t) G C(t), then 
by the results in Section 6.1, exit [ (I12|) has no required fundamental 
matrix] . 

(b) Assume that a = a(x)"/3(t) for some a{x) G C{x) and /3(t) G C(t). 
By the algorithms in [H [5], compute an irreducible matrix A such 
that A = (r{G)^G-^ for some G G GL„(fco). If ordoo(i) / 0, then 
by Proposition 13.261 exit [ ()12p has no required fundamental matrix]. 
Otherwise, let ^0 = ^U=oo and f3i{t),--- , f3n{t) be the eigenvalues 
of io. 
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(c) Goto Step (di) if the (3i{t) are conjugate and goto Step (^2) if = 
••• = /3„,(t) G C{t). In other cases, by Lemma 13.221 and Proposi- 
tion [3?26l exit [ (jl2p has no required fundamental matrix]. 

(di) If (Ai) yields no rational solutions, then exit [ (|12|) has no required 
fundamental matrix]. Otherwise, suppose that we find G G GL„(A;) 
such that 

a(G)a(rE)diag(/3i(t),--- ,pn{t)) = AG. 

Then B := G^^BG — G^^6{G) is of diagonal form. Compute a funda- 
mental matrix H of 

a{Y) = a(x)diag(/5i(t), • • • ,Pn{t))Y, 5{Y) = BY. 

Return [GH is a required fundamental matrix of ()12p ]. 

(^2) If we can compute a matrix G G GL„,(A;o) such that a{G)a{x)l3i{t) = 
AG then let 

B = G-'BG - G-'5{G) - ^-^^^xln G QUC{t)), 

Pi [t) 

else exit [ (jl2p has no required fundamental matrix]. If we can find a 
fundamental matrix H of d{Y) = BY whose entries are hyperexpo- 
nential over C{t), then reurn [GHhf3\{tY is a required fundamental 
matrix of ()12|) ] where h satisfies cr(/i) = a{x)h and = 0. Other- 
wise, exit [(|12p has no required fundamental matrix]. 

Decision Procedure 2 Compute a fundamental matrix of (jl3p whose en- 
tries are the interlacing of hypergeometric vectors over Uq if it exists. 

(a) If det(A) / (-l)"-i^a(x)/5(t) holds for any 5 G fc, /3(t) G C(t) 
and Q!(x) G C(x) that is standard with respect to cr", then exit [ ()13p 
has no required fundamental matrix]. 

(6) Expand det(74) as a series at a; = 00 : 

det(A) = + + • • • 

where P(t),(3i(t) G C(i) and m G Z. Suppose that x = p is not a 
pole of the entries of and that det(^o) / where = ;gpyU=p- 

Use (^42) to find all hypergeometric solutions of C7{Z) = AQ{nx)Z. By 
interlacing, we get all solutions of a^{Y) = AqY of the form Wihi. 
Denote these soluions by Wihi, • • • , Wdhd where Wi G C(x)"' and hi 
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satisfies cr'^{hi) = aihi for some ai G C(x). If there is io G {!,..., d} 

h A 

such that cj"(y) = ^„ ^^^^^ y has a rational solution in /cq, then 
let A(x) = ^^^17^, else exit [(USD has no required fundamental ma- 

trix]. Let jo be the least integer such that a^iY) = X(S+^^Xa(iy^ ^ 
rational solution in k^. If we can compute G G GL„(A;o) such that 

(7(y)/?(t)diag(A(x + jo),-- - ,A(x+jo + ?^- 1)) =^G, 

then let B = G^^BG - G''^6{G). So B is of diagonal form and by the 
same process as in Step (di) of Decision Procedure 1, we can compute 
a required fundamental matrix of (I13p . Otherwise, by the results in 
Section [3.3.21 exit [ (jl3p has no required fundamental matrix]. 

We can decide whether (|12p has liouvillian solutions or not as follows. If 
we can compute hypergeometric solutions over k of (jl2p by Decision Proce- 
dure 1, then we are done. Otherwise, consider the system (jl3p . If we can 
compute liouvillian solutions over /cq of (fT3j) by Decision Procedure 2, then 
by the results in Section 12.31 we can compute liouvillian solutions over ko 
of ()12p and we are done. Otherwise (I12p has no liouvillian solutions. 
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